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RECENT  METHODOLOGICAL  DEVELOPMENTS  IN  MAGNITUDE 
DETERMINATION  AND  YIELD  ESTIMATION  WITH  APPLICATIONS 
TO  SEMIPALATINSK  EXPLOSIONS 


R.-S.  Jih  and  R.  A.  Wagner 
Teledyne  Geotech  Alexandria  Laboratory 
314  Montgomery  Street 
Alexandria.  VA  22314-1581 

SUMMARY 

This  report  includes  two  parts.  The  first  section  discusses  in  detail  the  research 
that  was  done  since  the  submittal  of  our  first  annual  report  (GL-TR-90-0107),  and  the 
second  part  gives  a  perspective  overview  of  the  whole  project.  Note  that  the  same 
contract  number  (F19628-89-C-0063)  is  also  used  by  two  other  totally  independent 
tasks.  This  volume  covers  only  the  work  performed  under  Task  1. 

An  improved  magnitude  determination  procedure  is  described  and  tested  in  Sec¬ 
tion  I.  This  procedure  accounts  for  the  near-source  focusing/defocusing  effects  as 
well  as  the  receiver  effects  with  empirically  determined  correction  terms.  Although  no 
a  priori  geophysical  information  is  required  in  deriving  these  receiver  and  near-source 
terms,  the  inferred  corrections  turn  out  to  show  fair  correlation  with  the  tectonics 
underneath  the  receivers  as  well  as  the  visible  geological  structures  near  the  source 
region.  For  79  out  of  82  Semipalatinsk  events  in  our  WWSSN  database,  the  new 
scheme  provides  more  stable  measurements  across  the  whole  recording  network 
with  a  reduction  in  the  fluctuational  variation  by  a  factor  of  up  to  3.  The  3  events 
which  do  not  show  significant  improvement  could  have  been  detonated  in  environ¬ 
ments  with  different  focusing  patterns.  The  scatter  in  the  network- averaged  based 
on  the  new  scheme  versus  log(yield)  is  smaller  than  that  for  conventional  GLM  or 
LSMF  if  the  standard  and/or  the  rounding  errors  in  the  and  Soviet-published 
yields  (Bocharov  et  at. ,  1989)  are  included  in  the  regressions. 


Section  II  summarizes  the  research  accomplished  during  the  contract  period, 
including  those  major  results  already  reported  in  our  first  annual  report  (GL-TR-90- 
0107)  as  well  as  those  discussed  in  Section  I.  The  most  important  developments  in 
the  seismic  yield  estimation  methodology  under  this  project  are; 

•  a  refined  determination  scheme  (cf.  Section  I), 

•  an  algorithm  “MLE-CY”  which  utilizes  the  bounded  yields  in  the  -yield  regres¬ 
sion  based  on  the  maximum-likelihood  approach, 

•  an  algorithm  “DWLSQ”  which  permits  both  variables  to  be  imprecise  due  to  either 
rounding  or  standard  measurement  errors. 

We  have  fully  tested  these  techniques  with  Soviet-published  yields  and  the  WWSSN 
database  established  at  Geotech,  and  the  results  all  appear  to  be  very  encourag¬ 
ing  in  improving  our  remote  monitoring  capability. 

Also  included  in  this  report  is  a  complete  listing  of  192  event  m^,  values  measured 
off  21547  WWSSN  recordings  (Appendix  A).  We  have  updated  the  “yield- 
dependent"  test  site  bias  estimate  using  these  GLM91 A  values  (Appendix  B). 
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SECTION  I 


A  REFINED  NETWORK  DETERMINATION  SCHEME 
INCORPORATING  NEAR-SOURCE  EFFECTS 


Rong-Song  Jih  and  Robert  A.  Wagner 
Teledyne  Geotech  Alexandria  Laboratories 
314  Montgomery  Street 
Alexandria.  VA  22314-1581 


1.1  ABSTRACT 

An  improved  magnitude  determination  procedure  is  presented  to  account  for  the 
empirical  near-source  focusing/defocusing  effects.  For  79  out  of  82  Semipalatinsk 
events  in  our  WWSSN  database,  the  new  scheme  provides  more  stable  measure¬ 
ments  across  the  whole  recording  network  with  a  reduction  in  the  fluctuational  varia¬ 
tion  by  a  factor  of  up  to  3.  The  standard  error  in  the  resulting  network-averaged  is 
typically  around  0.02  m.u.,  similar  to  that  of  RMS  Lg  based  on  in-country  regional 
recordings.  The  3  events  which  do  not  show  significant  improvement  could  have  been 
detonated  in  environments  with  different  focusing  patterns.  The  scatter  in  the 
network-averaged  based  on  the  new  scheme  versus  log(yield)  is  smaller  than  that 
for  conventional  GLM  or  LSMF  if  the  standard  and/or  the  rounding  errors  in  the 
m^  and  Soviet-published  yields  are  included  in  the  regressions. 


I.2  INTRODUCTION 

The  main  problem  with  the  conventional  is  that  it  is  a  rather  nebulous  parame¬ 
ter;  simply,  it  is  a  function  of  the  largest  peak-to-peak  amplitude  in  the  first  few 
seconds  of  P  wave  motion  with  adjustment  for  the  period  of  the  arriving  phase.  The 
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parameter  was  adapted  from  the  need  to  order  systematically  the  size  of  earth¬ 
quakes.  The  measure  itself  has  inherent  impreciseness  as  the  measure  is  not  related 
to  the  physics  of  the  source  per  se  but  is  the  largest  constructive  interference  of 
waves  originating  at  the  source,  source  region,  propagation  path,  receiver  region,  and 
receivers  (Butler,  1981;  Johnson,  1981).  To  relate  the  m^,  to  the  seismic  yield,  all 
effects  not  due  to  the  source  must  naturally  be  corrected.  It  is  often  difficult,  however, 
to  separate  these  effects.  In  fact,  the  effects  of  source  and  propagation  are  often 
indistinguishable,  and  unless  one  is  known  the  other  can  not  be  uniquely  determined 
(Johnson,  1981).  Consequently,  it  was  reported  to  be  difficult  to  make  measure¬ 
ments  that  are  internally  consistent  within  0.1  with  the  conventional  (Bache, 
1982),  simply  because  some  of  the  aforementioned  effects  are  not  accounted  for 
accurately. 

Von  Seggern  (1973)  showed  that  including  station  corrections  typically  halves  the 
standard  deviation  of  from  North  American  LRSM  stations  recording  NTS  events. 
Even  better  results,  with  the  standard  deviation  reduced  by  a  factor  of  3  or  so,  can  be 
obtained  for  a  network  with  all  stations  beyond  30°.  Applying  station  corrections  to  the 
determination  or  the  network  spectra  averaging  has  become  a  standard  procedure 
in  this  community  (e.g. ,  Lilwall  etal..  1988;  Murphy  etai,  1989;  Sykes  and  Ekstrom, 
1989;  Jih  and  Shumway,  1989;  and  many  others).  The  station  effects  are  strongly 
dependent  on  azimuth  (Chang  and  von  Seggern,  1980),  which  led  Bache  (1982)  and 
many  others  to  believe  that  statistical  station  corrections  will  not  be  nearly  so  effective 
in  reducing  variance  in  the  multiple  source  region  problem, 

Marshall  et  al.  (1979)  attempted  to  correct  several  important  factors  that  can  bias 
m^,.  They  used  bulletin  log(A/T)  data.  These  data  are  corrected  for  receiver-station 
attenuation  differences,  and  the  resulting  magnitude  is  called  m2.  Correcting  m2  for 
source-region  attenuation  gives  m3,  and  correcting  m3  for  source  depth  gives  rriQ. 
(The  network  averages  of  these  m^s  are  denoted  by  ,  m2,  m3,  and  oIq,  respec¬ 
tively.)  The  major  change  in  this  scheme  is  associated  with  the  source-region 
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correction,  which  can  be  as  different  as  0.4  m.u.  between  sites  or  a  factor  of  about 
2.5  in  yield  estimates.  Essentially  this  approach  is  based  on  the  the  discovery  by 
Marshall  and  Springer  (1976)  and  Douglas  et  al.  (1981)  that  LRSM  amplitude  residu¬ 
als  correlate  with  P„  velocity  near  the  stations,  and  under  the  assumption  that  such 
correlation  between  the  attenuation  and  P„  is  valid  elsewhere  as  well.  However,  it 
turns  out  that  the  standard  deviation  of  the  niQ  is  not  less  than  that  for  the  m,  and  m2 
from  the  same  data  set,  indicating  that  the  attenuation  correction  in  Marshall  et  al. 
(1979)  is  not  correlated  with  the  residuals  from  these  stations  (Bache,  1982). 

It  is  obvious  that  the  only  way  to  reduce  the  statistical  fluctuation  is  to  obtain  fun¬ 
damental  causal  knowledge  of  the  focusing  and  defocusing  beneath  the  source  and 
receiver.  We  expect  teleseismic  P-wave  amplitudes  to  vary  as  source  location 
changes  within  a  test  site.  The  residuals  (with  respect  to  the  best-fitting  -yield 
curve)  of  NTS  events  show  systematic  trends  that  are  consistent  with  local  tectonics 
(Minster  et  al. ,  1981).  At  Yucca  Flat,  the  residuals  are  positive  to  the  west  and  nega¬ 
tive  to  the  east  of  the  north-south  trending  normal  fault  system  that  bisects  the  valley. 
At  Pahute  Mesa,  the  spatial  pattern  is  less  clear,  but  the  residuals  tend  to  be  negative 
toward  the  center  of  the  buried  Silent  Canyon  Caldera  and  positive  toward  the  edges. 
An  attractive  explanation  is  that  these  variations  are  due  to  focusing/defocusing  effects 
that  are  not  averaged  out  over  the  network,  although  the  possibility  of  systematic 
source-coupling  difference  has  not  been  eliminated. 

Figure  1  illustrates  the  residual  patterns  of  E.  Kazakh  explosions  as  seen  from 
various  directions.  For  each  event,  the  residual  is  the  (maximum-likelihood)  average 
of  all  station  residuals  {viz,  with  network  and  the  station  term  removed)  in  that 
quadrant.  All  three  subsites  exhibit  very  different  azimuthal  variation.  For  instance. 
Murzhik  events  are  enhanced  in  the  NE  and  SW  directions  and  reduced  in  the  NW 
and  SE  directions,  whereas  Degelen  events  are  reduced  towards  the  SW  direction. 
There  seems  to  be  some  weak  distinction  between  NE  and  SW  subregions  of  Balapan 
test  site  along  SE  and  NW  directions.  Figure  2  gives  the  azimuthal  pattern  with  the 
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four  quadrants  rotated  by  45°.  The  initial  P  waves  from  the  three  adjacent  test  sites 
have  virtually  the  same  incident  angle  at  any  particular  teleseismic  station,  and  any¬ 
thing  in  common  across  all  events  (such  as  the  crustal  amplification  as  well  as  the 
upper  mantle  attenuation  underneath  the  receiver)  would  have  been  lumped  into  the 
constant  station  term.  Thus  the  station  residuals  averaged  over  all  events  from  the 
same  test  site  would  correlate  very  little  with  the  receiver.  Instead,  they  should  reveal 
more  site-dependent  information  about  the  focusing/defocusing  pattern  underneath  E. 
Kazakhstan. 

In  this  study,  we  present  an  improved  scheme  to  determine  the  network  with 
both  the  station  terms  and  near-source  focusing/defocusing  effects  corrected.  Exam¬ 
ples  are  given  to  illustrate  that  such  procedure  can  reduce  the  random  fluctuation  in 
the  station  m*  values  to  about  0.15  m.u.  or  lower.  It  is  shown  that,  by  applying  this 
scheme  to  worldwide  explosions,  it  is  possible  to  have  a  consistent  base  line  in 
estimating  the  absolute  magnitudes,  which  is  crucial  in  estimating  the  test  site  bias, 
while  the  precision  in  the  resulting  network  values  can  be  maintained  as  well  as 
could  be  achieved  by  the  single-test-site  approach. 
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Figure  1.  Tfio  m.,  r-jslduals  of  F.  Kazakh  explosions  as  seen  from  various  directions  Tor  each 
event,  ttie  residual  is  the  (maximum-likelihood)  average  of  all  slaiion  residuals  (vi;!  with  network 
and  Ihe  slahon  term  removed)  in  that  quadrant.  All  three  subsites  exhibit  very  dilfcrenl  azimuthal 
variation  I'or  instance',  Murzhik  events  are  enhanced  in  the  NE  and  SW  diroclions,  and  reduced  in 
Ihe  NW  and  SE  directions,  whereas  Degelen  events  are  reduced  towards  the  SW  direction  There 
seems  to  be  some  weak  distinction  between  NE  and  SW  subregions  ol  Balapan  lest  site  along  SE 
and  NW  directions 
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1.3  DEFINITION  OF  A  NEW  MAGNITUDE,  m2  g 

The  conventional  definition  of  station  magnitude  is  computed  as 

=  logio(A/T)  +  B(A)  [1] 

where  A  is  the  displacement  amplitude  (in  nm)  and  T  is  the  predominant  period  (in 
sec)  of  the  P  wave.  The  B(A)  is  the  distance-correction  term  that  compensates  for  the 
change  of  P-vjawe  amplitudes  with  distance  {e.g.,  Gutenberg  and  Richter,  1956;  Veith 
and  Clawson,  1972).  in  [1]  is  also  denoted  as  in  Marshall  et  al.  (1979).  The 
ISC  bulletin  is  just  the  network  average  of  these  raw  station  rrii^  values  without  any 
further  adjustment. 

Consider  explosions  detonated  at  Np  source  regions  that  are  recorded  at 
some  or  all  of  Ng  stations.  The  GLM91A  network  mf,  (cf.  Appendix  A)  is  the 
(maximum-likelihood)  average  of  the  "station-corrected”  magnitudes: 

m2.2(i.j)  =  m,  (i,j)  -  S(j)  [2] 

where  S(j)  is  the  “statistical"  receiver  correction  at  the  j-th  station.  In  Marshall  et  al. 
(1979),  a  priori  information  about  the  velocity  underneath  each  station  is  used  to 
determine  its  associated  “deterministic”  receiver  correction,  S(j),  and  the  resulting 
magnitude  is  called  m2.  Our  receiver  corrections  (Figure  3),  however,  are  inferred 
jointly  from  a  suite  of  event-station  pairs,  and  no  a  priori  geophysical  or  geological 
condition  is  assumed  (and  hence  the  different  notation  m2  2)-  It  turns  out  when  the 
azimuthal  coverage  is  broad  enough,  receiver  corrections  derived  by  such  statistical 
approach  do  reveal  the  average  tectonic  structure  underneath  the  recording  stations, 
as  many  earlier  studies  have  reported  (e.g.,  North,  1977):  the  station  terms  are  posi¬ 
tive  in  shields  regions  such  as  Australia,  India,  Canada,  and  Scandinavia:  and  they  are 
negative  in  the  east  Africa  rift  valleys,  island  arcs  (e.g. ,  Japan  and  Taiwan),  and 
Indonesia.  The  high  correlation  between  the  tectonic  type  and  the  station  terms  sug¬ 
gests  that  the  station  corrections  do  reflect  the  upper  mantle  conditions  underneath  the 
receivers.  The  result  also  supports  the  claim  of  a  marginal  superiority  of  WWSSN 
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over  ISC  data.  For  instance,  Pacific  island  arcs  are  believed  to  have  high  attenuation, 
low  P„  velocities  as  well  as  negative  station  terms.  However,  North’s  (1977)  station 
corrections  based  on  38316  ISC  recordings  of  earthquakes  around  the  world  do  not 
show  this  phenomenon  because  ISC  does  not  have  so  wide  an  azimuthal  coverage 
and  uniform  spatial  sampling  as  does  WWSSN. 
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WWSSN  &  ISC  STATION  CORRECTIONS 
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Figure  3.  The  station  terms  derived  with  WWSSN  and  ISC  recordings.  Our  station  terms  (top; 
used  in  this  study)  are  based  on  the  GLM/MLE  joint  inversion  o(  192  worldwide  explosions  recorded 
at  122  “good"  WWSSN  stations.  Only  paths  within  20  and  95  degrees  are  used  l  or  each  station, 
our  GLM/LSMF  joint  inversion  scheme  puts  any  component  constant  across  all  events  as  the  "sta¬ 
tion  term",  similar  to  the  Douglas'  (1966)  LSMF  approach.  Both  GLM  (top)  and  LSMF  (middle)  sta¬ 
tion  corrections  exhibit  a  good  correlation  with  the  tectonics.  The  high  correlation  between  the  tec¬ 
tonic  type  and  the  station  terms  suggests  that  the  station  corrections  do  relied  the  upper  mantle 
conditions  underneath  the  receivers.  The  result  also  supports  the  claim  of  a  marginal  superiority  of 
WWSSN  over  ISC  data,  such  as  the  38316  ISC  recordings  used  by  North  (1977)  (bottom) 
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We  now  define  a  new  magnitude.  012^,  to  account  for  the  near-source  focusing 
and  defocusing  effects: 

m2.9(i.j)  =  ^”1  (i.j)  -  S0  -  F(k(i).j)  =  /t^.2  O.j)  -  F(k(i).j)  [3] 

At  the  j-th  station,  F(k(*).j)  is  a  constant  for  all  events  detonated  in  the  same  k-th  "geo¬ 
logically  and  geophysically  uniform  region”.  Partitioning  a  single  nuclear  test  site  into 
several  “regions”  may  be  necessary  in  order  to  account  accurately  for  the 
focusing/defocusing  effects.  This  m2g  is  very  similar  to  the  m3  in  Marshall  etal. 
(1979)  except  that,  again,  a  priori  attenuation  information  of  the  source  region  is  used 
in  Marshall  et  al.  (1979)  to  determine  the  correction  term,  whereas  we  invert  for  the 
near-source  effects  from  the  data  empirically.  (The  correlation  between  our  statistical 
focusing/defocusing  corrections  and  the  geological  structures  will  be  verified  in  a  later 
section.)  As  a  result,  the  source-region  corrections  used  by  Marshall  etal.  (1979)  are 
constants  (for  all  explosions  in  the  same  source  region)  regardless  of  the  location  of 
the  seismic  stations,  whereas  our  near-source  corrections  are  dependent  on  the 
source-station  paths. 

Table  1  lists  the  station  corrections  (which  are  invariant  for  any  explosion  from 
any  test  site  at  any  direction)  as  welt  as  the  near-source  corrections  associated  with 
each  subsite  of  Soviet’s  Semipalatinsk  nuclear  test  ground.'  Figure  4  plots  out  the 
lower-hemisphere  equal-area  projection  of  these  "secondary  corrections".  The  spatial 
map  of  these  corrections  are  shown  in  Figure  5. 


'Tho  noar-soufco  corroctioos  tor  oltier  tost  silos  can  be  stniilarly  derivod  Our  explosion  d.il.o  sol  needs  to  be  expandrjd. 
however,  to  warrant  a  reasonable  partitioning  at  other  source  regions 


EQUAL-AREA  PLOTS  OF  STATION  MEAN  mb(Pmax)  ANOMALIES 
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Figure  4.  Lower  heniisphore  equal-area  projection  o!  slalion-correcled  anomalies  (or 

Soviet’s  3  test  sites  in  Eastern  Kazakhstan.  For  each  station,  the  mean  m*  anomaly  is  dolmed  as 
the  mean  residual  averaged  over  all  explosions  from  the  same  test  site.  The  diflerent  patterns  more 
or  less  reflect  (he  focusing/defocusing  underneath  the  source  region. 
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MEAN  STATION-CORRECTED  mb{Pmax)  ANOMALIES 
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Figure  5.  Same  as  Figure  4  except  the  mean  mb  anomalies  are  superim(X)secl  by  the  map  We 
propose  to  regard  these  (site-dependent)  mean  station  mb  anomalies  as  the  secondary  correc 
lions"  to  be  applied  in  the  procedure  of  magnitude  determination 
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Table  1.  Receiver  and  Near-; 


Station  Term 


Near-source  Term.  F 


Code 


Balapan 


Degelen 


Murzhik 


AAE 


-0.352 


-0.387 


-0.076 


-0.127 


AAM 


0.205 


0.202 


-0.073 


-0.247 


ADE 


-0.044 


AFI 


-0.135 


AKU 


-0.048 


0.300 


0.311 


0.212 


ALQ 


-0.028 


ANP 


-0.349 


-0.167 


0.402 


0.183 


ANT 


0.040 


AQU 


-0.133 


-0.195 


0.039 


-0.024 


ARE 


0.223 


ATL 


0.141 


ATU 


0.128 


0.191 


-0.156 


0.034 


BAG 


-0.027 


0.076 


-0.009 


-0.103 


BDF 


0.067 


BEC 


-0.114 


0.090 


0.058 


-0.092 


BMP 


-0.219 


BKS 


0.089 


-0.014 


0.101 


0.229 


BLA 


0.057 


-0.233 


-0.177 


-0.299 


BOG 


0.032 


BOZ 


0.188 


-0.325 


-0.025 


-0.180 


Corrections  ol  WWSSN  Stations 


Station 


Longitude 

Latitude 

Description 

38.766 

9.029 

Addis  Ababa.  Ethiopia 

-83.656 

42.300 

Ann  Arbor,  Michigan 

138.709 

-34.967 

Adelaide,  south  Australia 

-171.777 

-13.909 

Afiamalu,  Samoa  Islands 

-18.107 

65.687 

Akureyri,  Iceland 

-106.457 

34.943 

Albuquerque,  New  Mexico 

121.517 

25.183 

Anpu,  Taiwan 

-70.415 

-23.705 

Antofagasta,  northern  Chile 

13.403 

42.354 

Aquila,  central  Italy 

-71 .491 

-16.462 

Arequipa,  southern  Peru 

-84.338 

33.433 

Atlanta,  Georgia 

23.717 

37.972 

Athens  Univ.,  Greece 

120.580 

16.411 

Baguio  City,  Luzon  Island 

-47.903 

-15.664 

Brasilia  Array,  Brazil 

-64.681 

32.379 

Bermuda-Columbia,  Atlantic 

-79.558 

8.961 

Balboa  Heights,  Panama 

-122.235 

37.877 

Byerly,  central  California 

-80.421 

37.211 

Blacksburg,  West  Virginia 

-74.065 

4.623 

Bogota,  Colombia 

-111.633 

45.600 

Bozeman.  Montana 

Table  1 .  Receiver  and  Near-source  Corrections  of  WWSSN  Stations 


Station  Term 


Near-source  Term,  F 


Code 

S 

BUL 

0.003 

CAR 

0.190 

CHG 

-0.140 

CMC 

-0.178 

COL 

0.065 

COP 

0.127 

COR 

0.155 

CTA 

0.153 

DAG 

0.036 

DAL 

0.202 

DAV 

-0,320 

DUG 

0.149 

EIL 

0.004 

EPT 

-0.023 

Longitude 

Latitude 

28.613 

-20.143 

-66.928 

10.507 

98.977 

18.790 

-115.083 

67.833 

-147.793 

64.900 

12.433 

55.683 

-123.303 

44.586 

146.254 

-20.088 

-18.770 

76.770 

-96.784 

32.846 

125.575 

7.088 

-112.813 

40.195 

34.950 

29.550 

-106.506 

31.772 

-3.205 

55.317 

-90.370 

38.802 

-90.426 

37.984 

-53.533 

69.250 

Description 


-0.147 


-0.003 

-0.175 


0.240 


0.114 


0.181 


-0.003 


0.132 


-0.072 


-0.052 


0.106 


0.375 


0.188 


0.159 


0.183 


0.003 


0.086 


0.045 


0.602 


0.051 


-0.276 


0.172 


-0.073 


-0.264 

-0.053 

0.038 

0.371 

0,352 

-0.117 

-0.228 

-0.103 

-0.042 

0.162 

-0  294 

-0.093 

-0.038 

0.045 

0.094 

0.015 

0.038 

-0.016 

r  ‘ 

0.336 


-0.051 


-77.067 

-90.300 


38.900 

-0.733 


Bulawayo,  Rhodesia 


Caracas,  Venezuela 


Chiengmai,  southeast  Asia 


Copper  Mine,  Canada 


College  Outpost,  Alaska 


Copenhagen,  Denmark 


Corvallis,  Oregon 


Charters  Towers,  Australia 


Danmatkshavn,  Greenland 


Dallas,  central  Texas 


Davao,  Mindanao  Island 


Dugway,  Utah 


Eilat,  Arabic  Peninsula 


El  Paso,  Texas-Mexico  border 


Eskdalemuir,  Scotland 
Florissanl.  eastern  Missouri 


French  Village,  eastern  Missouri 


Godhavn,  western  Greenland 


Georgetown,  Washington  D  C, 
Galapagos  Islands 
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Table  1.  Receiver  and  Near-source  Corrections  of  WWSSN  Stations 


Station  Term 


Near-source  Term,  F 


Longitude 

Latitude 

Description 

-105.371 

39.700 

Golden,  Colorado 

26.573 

-33.313 

Grahamstown,  southern  Africa 

-116.805 

35.302 

Goldstone,  central  California 

144.912 

13.538 

Guam,  Mariana  Islands 

114.172 

22.304 

Hong  Kong 

31.342 

29.858 

Helwan,  Arabic  Peninsula 

-67.986 

46.162 

Houlton,  New  Brunswick 

159.947 

-9.432 

Honiara,  Solomon  Islands 

28.996 

41.046 

Istanbul,  Turkey 

-99.802 

30.479 

Junction  City,  central  Texas 

35.197 

31.772 

Jerusalem,  Dead  Sea  region 

69.043 

34.541 

Kabul,  Afghanistan 

11.924 

78.918 

Kingsbay,  Svalbard  region 

27.007 

69.755 

Kevoa,  Finland 

-158.015 

21.423 

Kipapa,  Hawaii 

77.467 

10.233 

Kodaikanal,  India 

9.598 

59.649 

Kongsberg,  southern  Norway 

30.062 

69.724 

Kirkenes,  Sandinavia 

-21  983 

70.417 

Kap  Tobin,  eastern  Greenland 

74.333 

31.550 

Lahore,  India-Pakistan  border 

GOL 


GRM 


GSC 


GUA 


HKC 


HLW 


HN-ME 


0.114 


-0.075 


-0.072 


-0.059 


-0.130 


-0.096 


-0.204 


-0.285 


0.057 


-0.088 


-0.188 


-0.150 


0.175 


0  238 


0.148 


0.133 


0.011 


0.023 


-0.213 


-0.119 


0.196 


I 

0.046 

-0.004 

-0.208 

0.469 
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Table  1.  Receiver  and  Near-source  Corrections  of  WWSSN  Stations 


Station  Term 

Code 

S 

LEM 

-0.499 

LON 

-0.044 

LOR 

0.095 

LPA 

0.357 

LPB 

0.CS4 

LPS 

-0.068 

LUB 

0.191 

MAL 

-0.010 

MAN 

0.276 

MAT 

-0.188 

MDS 

-0.043 

MNN 

0.153 

MSH 

0.202 

MSO 

-0.045 

MUN 

0.172 

NAI 

-0.112 

NAT 

0.118 

NDI 

0.158 

NHA 

-0.134 

NIL 

-0.030 

Near-source  Term.  F 


Station 


0.123 

-0.535 

-0.119 

0.159 

0.077 

-0.362 

-0.165 

0.056 

-0.019 

-0.047  1 

— 

0.184 

0.076 

0026 

0.075 

-0.108 

-0.075 

Longitude 

Latitude 

Description 

107.617 

-6.833 

Lembang,  Java 

-121.810 

46.750 

Longmire,  Washington 

3.851 

47.267 

Lormes.  France 

-57.932 

-34.909 

La  Plata,  Uruguay 

-68.098 

-16.533 

La  Paz,  Peru-Bolivia  border 

-89.162 

14.292 

La  Palma,  Guatemala 

-101.867 

33.583 

Lubbock,  west  Texas 

-4.411 

36.728 

Malaga.  Straits  of  Gibraltar 

121.077 

14.662 

Manila.  Luzon  island 

138.207 

36.542 

Matsushiro,  Honshu,  Japan 

-89.760 

43.372 

Madison,  Wisconsin 

-93.190 

44.914 

Minneapolis,  Minnesota 

59.588 

36.311 

Meshed.  Iran-USSR  border 

-113.941 

46.829 

Missoula.  Mo.nlana 

116.208 

-31.978 

Mundaring,  western  Australia 

36.804 

-1.274 

Nairobi.  Kenya 

-35.033 

-5.117 

Natal,  Brazil 

77.217 

28.683 

New  Delhi,  northern  India 

109.212 

12.210 

Nhatrang,  southeast  Asia 

73.252 

33.650 

Nilore,  Pakistan 

i.  ..  . 
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Table  1.  Receiver  and  Near-source  Corrections  of  WWSSN  Stations 


Station  Term 


I  Code 

s 

NNA 

-0.171 

NOR 

-0.240 

NP  NT 

0  107 

NUR 

0.090 

OGD 

-0.167 

OXF 

0.293 

Near-source  Term,  F 


Murzhik  Longitude  Latitude 


Station 


Description 


0.473 


-0.061 


-0.076 


-0.214 


0.041 


-0.263 


-76.842  -11.988  Nana,  Peru 


-16.683  81.600  No rd,  north  coast  of  Greenland 


-119.372  76.252  North  Pole,  Queen  Elizabeth  Islands 

24.651  60.509  Nurmijarvi,  Finland 


-74.596  41 .088  Ogdensburg,  New  York 


-89.409  34.512  Oxford,  Mississippi 


-25.663  37.747  Ponta  Delgada,  Azores  Islands 


-70.685  -33.144  Peldehue,  Chile-Argentina  border 


PMG 

0.151 

0.101 

0.060 

0.253 

147.154 

-9.409 

Port  Moresby,  New  Guinea 

POO 

0.076 

-0.041 

0.082 

0.251 

73.850 

18.533 

Poona,  India 

PRE 

-0.089 

0.118 

-0.210 

-0.017 

28.190 

-25.753 

Pretoria,  south  Africa 

PTO 

-0.172 

-0.163 

-0.166 

-0.029 

-8.602 

41.139 

Porto,  Serro  Do  Portugal 

QUE 

-0  446 

0.484 

0.107 

0.190 

66.950 

30.188 

Quetta,  Pakistan 

QUI 


RAB 


R 


RCD 


RIV 


RK-ON 


0.023 


0.022 


.093 


0.370 


0.300 


-0.013 


-78.501  -0.200  Quito,  Ecuador 


152.170  -4.191  Rabaul,  New  Britain  region 


-159.773  -21.212  Rarotonga,  Cook  Islands  region 


-103.208  44.075  Rapid  city.  South  Dakota 


151.158  -33.829  Riverview,  SE  Australia 


-93.672  50.839  Red  Lake.  Ontario 


SCP 

-0  005 

0.037 

-0.129 

-0.229 

-77.865 

40.795 

Stale  College  Pennsylvania 

SDB 

0.053 

0.114 

0  130 

0.172 

13.572 

-14.926 

Sa  Da  Bandeira,  Angola 
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Table  1.  Receiver  and  Near-source  Corrections  of  WWSSN  Stations 


Station  Term 


Near-source  Term,  F 


Station 


Code  S  Balapan  Degelen  Mur2hiK  Longitude  Latitude 


Code 

S 

SEO 

-0.132 

SHA 

0.297 

SHI 

0.237 

SHK 

-0.286 

SHL 

-0.001 

SJG 

-0.129 

SNG 

-0.003 

SPA 

-0  630 

STU 

0.067 

,  TAB 

0.234 

TAU 

-0.137 

TOL 

0.169 

TRI 

-0,144 

TRN 

0.064 

TUC 

—J 

0.008 

Description 


0.144 

-0.266 


-0.027 


0.107 
-0.216 
-0  154 


-0.273 


-0.063 


0.106 


-0.027 


-0.387 


0.024 


-0.105 


-0.160 


-0.088 


-0.135 


0.027 


-0.143 

0.215 


-0.027 

0.016 


-0.014 

0.151 


-0.309 

0.009 


-0.141 

-0.055 


126.967  37.567  Seoul.  South  Korea 


-88.143  30.694  Spring  Hill,  Mississippi 


52.520  29.638  Shiraz,  southern  Iran 


132.678  34.532  Shiraki,  southern  Honshu,  Japan 


91.883  25.567  Shillong,  India-Bangladesh  border 


-66.150  18.112  San  Juan,  Puerto  Rico 


100.620  7.173  Songkhia,  Malay  Peninsula 


0.000  -90.000  South  Pole,  Antarctica 


9.195  48.772  Stuttgart,  Germany 


46.327  38.068  Tabriz,  Iran-USSR  border 


147.320  -42.910  Tasmania  Univ.,  Tasmania 


-4.049  39.881  Toledo.  Spain 

13.764  45.709  Trieste,  northern  Italy 


-61.403  10.649  Trinidad.  Trinidad 


-110.782  32.310  Tucson,  eastern  Arizona 


20  237  63  815  Umea,  Sweden 

-99.178  19.329  Nat  Univ  of  Central  Mexico 

-10.244  51.939  Valentia,  Eire 

174.768  -41  286  Wellington,  New  Zealand 

-71.322  42.385  Weston,  New  England 

17  100  -22  567  Windhoek,  South-West  Afnca 
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Table  2  lists  the  results  of  applying  the  receiver  and  near-source  corrections  (as 
shown  in  Table  1)  to  the  82  Semipalatinsk  explosions  in  our  database.  For  79  out  of 
82  Semipalatinsk  events  used  in  this  study,  the  final  a  is  typically  around  or  below  the 
same  level  that  a  "single-event  MLE  with  primary  correction  only”  could  achieve.  The 
only  three  events  which  do  not  show  reduction  in  the  variance  are  770730D.  730723B, 
and  88091 4B  (Table  2).  A  plausible  explanation  is  that  perhaps  these  three  events  are 
located  in  very  different  geological  or  geophysical  environments  from  other  events  in 
the  same  testing  area.  Or,  perhaps  the  26  WWSSN  stations  for  which  the  filmchips  of 
JVE  were  available  only  cover  a  small  portion  of  the  focal  sphere,  and  hence  the 
focusing/defocusing  effect  is  not  fully  accounted  for.  Most  of  the  /TT2  9  in  Table  2  have 
a  standard  error  around  0.02  m.u.,  which  is  about  the  same  as  that  for  RMS  Lg  values 
inferred  from  an  in-country  regional  network  {e.g.,  Israelson,  1991;  Hansen  etal., 
1990). 

Another  observation  is  that  the  resulting  values,  ffT2  9  are  essentially  the  same 
as  those  inferred  from  the  GLM  or  the  single-event  MLE,  /n2.2-  Figures  6  through  10 
plot  out  the  three  different  m^s  for  five  arbitrarily  selected  events.  The  solid  line  and 
the  dashed  lines  represent  the  mean  network-averaged  and  the  associated  stan¬ 
dard  deviation  of  station  m^.  The  upward  arrows  represent  the  lower  bounds  of  the 
station  which  came  from  those  clipped  measurements.  The  "Y”  symbols  are  the 
upper  bounds  of  the  station  m^j  which  are  resulted  from  those  noisy  measurements. 
Both  the  "uncensored"  (shown  in  filled  circles)  and  “censored”  station  mj,s  are  used 
in  computing  the  mean  station  residuals  with  the  maximum-likelihood  scheme 
described  in  Jih  and  Shumway  (1989). 
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Table  2.  Comparison  of  Network-Averaged  m.  with  Various  Corrections 


Event 


Date 


#  of  signals  Without  correction  Station  corrected 


661218M 


67091 6M 


670922M 


671122M 


690531 M 


691228M 


700721 M 


701104M 


711009M 


720902M 


51  2  1 


36  18  2 


35  20  1 


7  52  0 


30  21  0 


45  2  3 


38  12  1 


38  12  1 


38  6  2 


41  6  0 


27  9  3 


32  6  0 


29  10  2 


15  25  0 


HTi .  a 


5.747±0.038  0.281 


5.097±0.041  0.305 


5.03310.036  0.271 


4.16910.054  0.413 


4.96510.050  0.357 


5.66610.043  0.306 


5.18210.043  0.307 


5.26710.042  0.303 


5.32310.040  0.272 


5.29410.040  0.278 


5.16510.037  0.233 


5.34810.049  0.301 


5.16810.045  0.290 


4.61510.049  0.312 


^2.2- 


5.75310.033  0.239 


5.11010.036  0.268 


5.04810.029  0.214 


4.31810.039  0.299 


5.00610.039  0.276 


5.66510.035  0.250 


5.19910.033  0.236 


5.27910.032  0.232 


5.34110.031  0.210 


5.31110.030  0.208 


5.18710.028  0.174 


5.38310.036  0.224 


5.18610.030  0.193 


4.63510.038  0.243 


Near-source  corrected 


^2.9  >  ^ 


5.73810.022  0.161 


5.09510.018  0.137 


5.02910.017  0.125 


4.23110.013  0.099 


5.02610.015  0.110 


5.66010.018  0.125 


5.18410.018  0.125 


5.24910.020  0.145 


5.31910.015  0.099 


5.28710.013  0.086 


5.13610.016  0.100 


5.34110.021  0.127 


5.16310.017  0.111 


4.60210.017  0.107 


651121D 

48  12  1 

5.38410.035  0.271 

5.39410.024  0.188 

5.38110.019  0.152 

66021 3D 

51  2  10 

6.07310.038  0.305 

6.09210.027  0.218 

6.08810.014  0.114 

660320D 

49  6  8 

5.84810.041  0.322 

5.86510.030  0.239 

5.85310.009  0.074 

660507D 

9  23  1 

4.51710.043  0.250 

4.55910.031  0.181 

4.45610.014  0.078 

661019D 

51  8  5 

5.53410.035  0.276 

5.54210.025  0.201 

5.52510.014  0.111 

1 
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Table  2.  Comparison  of  Network-Averaged  with  Various  Corrections 


Event 


Date 


680929D 


690723D 


710425D 


711230D 


720328D 


72081 6D 


721210D 


770329D 


770730D 


780326D 


780422D 


780728D 


68061 9B 


691130B 


#  of  signals 

Without  correction 

Station  corrected 

Near-source  corrected 

,  a 

^2.2  •  ^ 

fTIz  g ,  a 

50  4  6 

5.610±0.035  0.275 

5.64210.026  0.202 

5.64110.018  0.138 

38  17  1 

5.172±0.041  0.306 

5.20110.029  0.220 

5.18610.013  0.100 

19  35  0 

4.533±0.053  0.392 

4.60510.040  0.295 

4.63410.025  0.186 

43  11  3 

5.49810.040  0.305 

5.52810.032  0.245 

5.53010.017  0.126 

37  3  0 

5.76410.052  0.327 

5.79310.042  0.267 

5.82610.020  0.127 

16  2  0 

5.52210.060  0.255 

5.55610.062  0.262 

5.55610.035  0.148 

28  15  0 

4.95510.048  0.313 

4.99710.034  0.226 

4.98410.021  0.141 

23  20  1 

4.90810.044  0.293 

4.93110.033  0.221 

4.92110.025  0.165 

30  6  5 

5.51210.045  0.287 

5.54310.033  0.209 

5.55510.029  0.183 

25  12  0 

4.97410.055  0.332 

5.01510.042  0.256 

4.99210.042  0.254 

21  14  0 

4.85810.051  0.301 

4.88110.049  0.287 

4.85510.052  0.306 

25  4  0 

5.46110.052  0.279 

5.51910.037  0.198 

5.47610.018  0.096 

21  7  0 

4.98510.057  0.301 

5.03210.044  0.231 

5.01010.023  0.120 

36  7  6 

5.50610.042  0.291 

5.52510.028  0.198 

5.49410.012  0.086 

36  20  1 

5.13810.033  0.250 

5.15610.025  0.193 

5.13110.012  0.089 

46  1  2 

5.89610.040  0.280 

5.89310.032  0.223 

5.86110.022  0.155 

28  3  2 

5.27610.047  0.270 

5.28510.035  0.203 

5.22910.026  0.151 

51  0  0 

5.95010.044  0.313 

5.97310.031  0.222 

5.94510.026  0.184 

31  18  1 

5.04510.043  0.301 

5.07510.035  0.247 

5.04310.027  0.192 

Table  2.  Comparison  of  Network-Averaged  with  Various  Corrections 


Event 


Date 


it  of  signals  Without  correction  Station  corrected 


m,.  o 


^2.2  ’  ^ 


Near-source  corrected 


1^2.9  ’  ® 


721102B 

42  0  15 

6.17310.045  0.339 

6.18310.034  0.256 

6.16010.023  0.177 

721 21 OB 

45  1  11 

6.00610.037  0.282 

6.01310.029  0.223 

5.98310.022  0.169 

730723B 

54  0  1 

6.17410.042  0.314 

6.20210.030  0.220 

6.17910.031  0.231 

731214B 

50  7  6 

5.75010.044  0.348 

5.76910.036  0.287 

5.74910.030  0.241 

750427B 

18  1  1 

5.45710.099  0.444 

5.48010.089  0.397 

5.43610.072  0.322 

760704B 

38  0  5 

5.81910.062  0.404 

5.84910.048  0.313 

5.82910.027  0.180 

761207B 

17  2  1 

5.57110.107  0.480 

5.61110.099  0.442 

5.55010.078  0.347 

78061  IB 

170  1 

5.88210.048  0.205 

5.87310.042  0.179 

5.80410.038  0.160 

78091 5B 

37  1  6 

5.82610.057  0.379 

5.85010.043  0.284 

5.82810.030  0.199 

790623B 

40  2  3 

6.01510.052  0.349 

6.07110.040  0.267 

6.06910.036  0.240 

790804B 

40  4  20 

6.06410.037  0.295 

6.08610.027  0.212 

6.10310.016  0.126 

1  791028B 

44  6  13 

5.90910.036  0.284 

5.94110.024  0.190 

5.93310.019  0.151 

791223B 

41  3  17 

6.11110.033  0.260 

6.12810.020  0.154 

6.13110.018  0.141 

80091 4B 

35  4  6 

6.00610.062  0.417 

6.04110.051  0.345 

6.05310.043  0.287 

811018B 

- 1 

41  3  7 

5.95410.039  0.279 

5.97610.029  0.210 

5.99610.023  0.164 

840526B  31  0  3  5.966±0.058  0.338 

880914B  25  0  1  6.004±0.037  0.191 

761 1 23B  22  0  0  5.577±0.075  0.354 

780829B  16  0  0  5.869)0  079  0.315 

781 1 29B  28  0  0  5.840±0.068  0.358 

790707B  30  0  0  5.734±0.048  0.261 


5.993±0.043  0.252 
6.032±0.023  0.117 
5.62610.065  0.305 
5.90510.075  0.302 
5.88010.054  0.287 
5.80010.043  0.236 


6.00910.022  0.130 
6.02610.033  0.168 
5.68710.053  0.249 
5.93610.045  0.182 
5.89510.026  0.138 
5.82710.031  0.169 
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Table  2.  Comparison  of  Network-Averaged  with  Various  Corrections 


Event 


Date 


#  of  signals  Without  correction  Station  corrected 


m. .  a 


072.2  •  ® 


Near-source  corrected 


OT2  9 ,  or 


79081 8B 

28  0  0 

6.02310.065  0.344 

6.08710.052  0.273 

6.08810.031  0.166 

791202B 

15  0  0 

5.80710.112  0.435 

5.87410.089  0.344 

5.89210.067  0.259 

801012B 

23  0  0 

5.80410.087  0.417 

5.82810.070  0.334 

5.87210.047  0.228 

801214B 

29  0  0 

5.87310.053  0.288 

5.91110.047  0.252 

5.93510.030  0.162 

801227B 

24  0  0 

5.85510.056  0.273 

5.89610.042  0.208 

5.89210.034  0.165 

810422B 

25  0  0 

5.82610.070  0.350 

5.86510.057  0.287 

5.92210.035  0.174 

810913B 

17  0  0 

6.01410.093  0.383 

6.02410.064  0.265 

6.06010.023  0.096 

811227B 

23  0  0 

6.18710.079  0.380 

6.20710.060  0.287 

6.18910.036  0.174 

820425B 

14  0  0 

5.92410.099  0.372 

5.94410.072  0.269 

5.98210.040  0.150 

820704B 

21  0  0 

6.07310.053  0.245 

6.08910.048  0.222 

6.09510.032  0.145 

821205B 

26  0  0 

6.04310.064  0.325 

6.09310.053  0.271 

6.10910.037  0.191 

83061 2B 

16  0  0 

5.92110.091  0.365 

5.94310.069  0,276 

5.93410.046  0.186 

831006B 

25  0  0 

5.88510.064  0.318 

5.94210.048  0.238 

5.93510.032  0.162 

831026B 

18  0  0 

5.93310.070  0.298 

5.94110.053  0.225 

5.99810.034  0.146 

840425B 

21  0  0 

5.85010.099  0.453 

5.89210.082  0.374 

5.90510.043  0.196 

84071 4B 

23  0  0 

5.92010.088  0.424 

5.99910.074  0.357 

6.05410.043  0.207 

841027B 

19  0  0 

6.14110.078  0.340 

6.15010.066  0.289 

6.19110.034  0.146 

841202B 

22  0  0 

5.63010.065  0.305 

5.69310.057  0.266 

5.72010.044  0.206 

841216B 

15  0  0 

5.91110.107  0.415 

5.99310.070  0.272 

6.04310.035  0.135 

841228B 

19  0  0 

5.85310.077  0.335 

5.91610.053  0.230 

5.94510.037  0.162 

85061 5B 


15  0  0 


6.01610.078  0.301  6.06910.049  0.191 


6.06010.035  0.134 


VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  710606K 


O  50  100  150  ZOO  Z50  300  350 


Ay.irnut/L  (degrees) 

Figure  6.  Scatter  plot  of  3  different  types  of  station  mj,s  for  Murzhik  explosion  710606  The  33 
good  recordings,  6  noise,  and  2  clips  are  shown  with  filled  circles,  Y-shaped  d:  .vn.vaio  '’'r.'.s 
upward  arrows,  respectively.  The  raw  station  mt,s  (top)  has  a  standard  devuiion  c  0  27  tnu 
Applying  the  “primary"  station  corrections  reduces  tlie  scatter  to  0  21  m  u  Applying  to-,-  proposed 
“secondary"  station  corrections  to  count  lor  the  near-source  focusing/detocusmg  would 

further  reduce  the  scatter  down  to  0.1  m.u.  The  dashed  lines  around  the  iietv^ork-a.eraged 
clearly  illustrate  the  remarkable  reduction  of  fluctuation  across  the  recordirrg  stations  Ihe  mean 
event  itself  is  not  significantly  changed,  however. 
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VAIilOUS  WWSSN  MAGNITUI)I‘:S  OF  KVKNT  710010K 


Figure  7.  Same  as  Figure  6  except  for  Murzhik  event  710619. 
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VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  711009K 


6.5 


Station  mj,'s  with  Mean  Receiver  Effects  Removed 


5.187  +/-  0.028  ,  a=  0.174  27.9.3 


6  0 


5.5  h- 


4.5  - 

^  Q  _ I _ 1 _ I - 1 - 1 - 1 - 1 - 

'  0  50  100  150  200  250  300  350 


Station  mUs  with  Receiver  &  Neai — source  Terms  Removed 

6.5  I - — - 

5.136  +/-  0.016  .  a-  0.100  27.9.3 

GO  - 

Y 

Y  Y 

5  ^  /  Y. 

5.0  ^  ^  —  •'*'  —  —  —r  J-  —  ___ 

4. 5 

4. o' - ' - ' - ' - ^ ^ ^ 

0  50  100  150  200  250  300  350 

Azimuth  (degrees) 


5.136  +/-  0.016  .  a-  0.100  27.9.3 


Y  Y 


Y 


Y 

Y. 


_ _ 


-  -  2^ 


J _ 1_ 


Figure  8.  Same  as  Figure  6  except  for  Murzhik  event  711009. 
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Aximuth.  (degrees) 


Figure  9.  Same  as  Figure  6  except  for  Degelen  event  660320,  The  near-source  correction  pro¬ 
posed  in  this  study  not  only  reduced  the  scatter  at  stations  that  reported  the  good  signals,  but 
also  improved  the  data  consistence  of  the  censored  recordings. 
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VARIOUS  WWSSN  MACNITUDES  OF  EVENT  710-1250 


Raxi)  Slaiion  mi's  (No  Co7~ractiori  AjipLiad) 


1.4  UNDERLYING  MODEL  AND  ADVANTAGES  OF  THE  NEW  MAGNITUDE 

We  now  examine  the  fundamental  difference  between  the  present  scheme  and 
the  previous  ones.  In  LSMF  and  the  standard  GLM  scheme  (Douglas,  1966;  Blandford 
and  Shumway,  1982;  Jih  and  Shumway,  1989;  Murphy  etal.,  1989),  it  is  assumed 
that  the  observed  station  mb(i,j)  is  the  sum  of  the  true  source  size  of  the  i-th  event, 
E(i),  the  receiver  term  of  the  j-th  station,  SQ),  and  the  random  noise,  v(i,j): 

"»b(i.j)  =  E(i)  +  S(j)  +  v(i,j)  [4] 

The  receiver  term,  S(j),  is  constant  with  respect  to  all  explosions  from  many  azimuths, 
and  hence  it  would  inherently  reflect  the  “averaged”  receiver  effect  —  provided  the 
paths  reaching  the  station  have  broad  azimuthal  coverage.  These  receiver  corrections 
correlate  with  the  upper  mantle  property  underneath  the  receivers  (North,  1977). 
When  world-wide  explosions  are  used,  the  standard  deviation  of  the  noise  v  in  [4]  is 
typically  about  0.3  m.u.  or  larger. 

If  LSMF  or  GLM  is  applied  to  events  within  a  smaller  area  of  source  region,  then 
the  a  could  reduce  to  0.15  or  0.2  m.u.  Unfortunately,  there  are  severe  drawbacks 
associated  with  such  “single-test-site  GLM”  approach.  First  of  all,  the  station  correc¬ 
tions  will  not  necessarily  represent  the  attenuation  underneath  the  receiver  side.  They 
could  be  contaminated  or  even  overwhelmed  by  the  near-source  effects  shared  by  the 
explosions  confined  in  a  narrow  azimuthal  range.  This  explains  the  phenomenon 
Butler  (1981)  and  Burdick  (1981)  reported  that  using  Soviet  explosions  exclusively 
may  fail  to  discern  the  attenuation  differential  between  the  eastern  and  western  U.S. 
Secondly,  when  the  “single-test-site  GLM"  inversion  is  applied  to  several  test  sites 
separately,  there  may  not  be  a  consistent  baseline  for  magnitude  comparison  or  abso¬ 
lute  yield  estimation,  since  the  station  terms  are  inherently  inconsistent. 

In  the  present  scheme  ([3]),  however,  we  reformulate  the  whole  model  as 

mf,(i,j)  =  E(i)  -(-  S(j)  +  F(k(i).j)  -t-  v(i,j)  [5] 

where  F(k(i),j)  is  the  correction  term  at  the  j-th  station  for  the  near-source 
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focusing/defocusing  effect,  which  is  constant  for  all  events  in  the  k-th  "geologically  and 
geophysically  uniform  region”.  For  each  seismic  station,  this  F  can  be  regarded  as  its 
azimuthal  variation  around  the  mean  station  term  S.  However,  as  we  already 
explained,  it  would  be  more  appropriate  to  consider  F  the  near-source  term  because 
the  back  azimuths  at  the  station  could  be  nearly  identical  for  adjacent  test  sites  (such 
as  Degelen  and  Murzhik),  and  yet  the  “F"  terms  could  be  very  different.  By  incorporat¬ 
ing  the  F  term  into  the  model,  the  o  for  world-wide  explosions  is  reduced  to  about  0.2, 
roughly  the  same  level  that  which  a  “single-test-site  GLM”  could  achieve.  Intuitively, 
the  present  scheme  (Equation  [5])  provides  a  more  detailed  (and  better)  model  than 
that  of  Equation  [4]  in  describing  the  whole  propagation  path  from  the  source  towards 
the  receiver.  Simply  put.  Equation  [4]  yields  a  stronger  fluctuation  in  the  source  terms, 
E,  as  well  as  a  larger  standard  deviation  of  v  because  each  term  in  the  right-hand  side 
of  Equation  [4]  would  have  to  "absorb"  part  of  the  missing  F  term  in  [5].  This  is 
exactly  the  same  reason  why  m22  has  smaller  variation  than  since  the  latter 
would  be  interfered  by  the  missing  station  term  S  in  Equation  [1]. 

For  actual  implementation,  the  present  scheme  can  be  replaced  with  an 
equivalent  multi-stage  procedure  as  follows.  First,  a  set  of  station  corrections  (the  so- 
called  "primary  corrections")  is  determined  with  one  GLM.  Then  the  "secondary 
correction"  at  each  station  is  defined  as  the  mean  of  all  residuals  of  all  events  from 
the  same  test  site  (or  the  same  geologic/geophysical  regime)  recorded  at  this  particu¬ 
lar  station. 


1.5  MAGNITUDE:YIELD  RELATIONSHIP  AT  SEMIPALATINSK  WITH 

To  further  demonstrate  that  rn2  g  would  provide  more  precise  yield  estimates,  the 
19  Semipalatinsk  explosions  for  which  the  yields  are  published  by  Bocharov  et  al. 
(1989)  are  used  as  a  test  case.  Table  3  gives  the  date,  various  values,  and  the 
associated  standard  errors.  Table  4  lists  the  Soviet-published  yields  and  the 


-  30- 


postulated  uncertainties.  We  assume  that  these  yields  are  subject  only  to  10%  stan¬ 
dard  errors  (S.E.)  and/or  the  rounding.  For  each  {mt,,  yield)  pair,  we  use  a  random 
number  generator  to  produce  a  perturbed  (m^,  yield)  pair  according  to  their  uncer¬ 
tainty  distribution.  A  standard  least-squared  regression  is  performed  for  each  data  set 
of  perturbed  samples.  The  procedure  is  repeated  for  several  hundred  iterations,  and 
the  resulting  calibration  curves  (i.e.,  the  straight  best-fitting  lines)  are  shown  as  the 
darkened  bundle  in  Figures  11  through  13.  The  detail  of  this  generalized  “doubly- 
weighted  least-squares  scheme"  is  discussed  in  Jih  (1991).  Here  we  only  summarize 
the  results  to  illustrate  the  advantages  of  /7T2  g  relative  to  the  more  conventional  source 
measures.  For  comparison,  regression  result  using  RMS  Lg  reported  at  NORSAR 
(Ringdal,  1990)  is  also  included  in  Tables  5  and  6  (Figure  14).  Note  that  for  the  Soviet 
JVE  shot  (88091 4B),  the  yield  is  assumed  to  be  119  kt  after  Gordan  (1988)  (see  also 
Sykes  and  Ekstrom,  1989;  Priestley  etai,  1990).  The  regression  result  based  on 
m29  has  a  smaller  scatter  around  the  mean  calibration  curve  (and  hence  a 
smaller  uncertainty  factor  in  the  yield  estimates)  than  those  based  on  and  r)72  2,  as 
expected.  It  has  a  precision  very  close  to  that  based  on  NORSAR  RMS  Lg  over  a 
wide  range  of  yields.  The  precision  is  also  very  similar  to  what  Patton  (1988)  found  for 
69  NTS  explosions  below  the  water  table  with  Lg  recorded  at  LLNL  digital  network. 
Note  that  the  uncertainty  factor  in  the  yield  estimates  is  yield  dependent.  It  is  smaller 
near  the  centroid  of  the  data  set  (i.e. ,  around  50  kt  in  our  case)  and  larger  at  both 
ends.  This  is  contrary  to  a  general  perception  that  we  might  know  yields  much  better 
at  higher  values  around  150  kt.  For  yields  below  10  kt  there  is  no  data  point  in 
NORSAR’s  RMS  Lg  data  set,  and  hence  there  is  a  much  larger  uncertainty  than  that 
based  on 
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Table  3.  Various  Network-Averaged  of  19  Special  Events 

Event 

#  of  signals 

Without  correction 

Station  corrected 

Near-source  corrected 

Date 

.  a 

^2  2  •  ® 

^2.9  <  ® 

651121D 

48  12  1 

5.384±0.035  0.271 

5.39410.024  0.188 

5.38110.019  0.152 

66021 3D 

51  2  10 

6.07310.038  0.305 

6.09210.027  0.218 

6.08810.014  0.114 

660320D 

49  6  8 

5.84810.041  0.322 

5.86510.030  0.239 

5.85310.009  0.074 

670922M 

35  20  1 

5.03310.036  0.271 

5.04810.029  0.214 

5.02910.017  0.125 

680929D 

50  4  6 

5.61010.035  0.275 

5.64210.026  0.202 

5.64110.018  0.138 

690723D 

38  17  1 

5.17210.041  0.306 

5.20110.029  0.220 

5.18610.013  0.100 

691130B 

51  0  0 

5.95010.044  0.313 

5.97310.031  0.222 

5.94510.026  0.184 

691228M 

45  2  3 

5.66610.043  0.306 

5.66510.035  0.250 

5.66010.018  0.125 

710425D 

37  3  0 

5.76410.052  0.327 

5.79310.042  0.267 

5.82610.020  0.127 

710606M 

38  6  2 

5.32310.040  0.272 

5.34110.031  0.210 

5.31910.015  0.099 

711009M 

27  9  3 

5.16510.037  0.233 

5.18710.028  0.174 

5.13610.016  0.100 

71 1021 M 

32  6  0 

5.34810.049  0.301 

5.38310.036  0.224 

5.34110.021  0.127 

72021  OB 

34  8  2 

5.29710.042  0.278 

5.32910.029  0.195 

5.28910.018  0.121 

72081 6D 

23  20  1 

4.90810.044  0.29? 

4.9?'' 10  03?  0.221 

4.92110.025  0.165 

720902M 

15  25  0 

4.61510.049  0.312 

4.63510.038  0.243 

4.60210.017  0.107 

721102B 

42  0  15 

6,17310,045  0,339 

6.18310.034  0.256 

6.16010.023  0.177 

721 21  OB 

45  1  11 

6.00610.037  0.282 

6.01310.029  0.223 

5.98310.022  0.169 

88091 4B 

25  0  1 

6.00410.037  0.191 

6.03210.023  0.117 

6.02610.033  0.168 
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Table  6.  95%  Confidence  Scatter  in 

(assuming  yields  are  subject  to  rounding  and  10%  S.E.) 

used 

t  „  — — . , 

1  kt 

10  kt 

50  kt 

100  kt 

150  kt 

m, 

0.27 

0.16 

0.14 

0.15 

0.17 

^2.2 

0.24 

0.14 

0.11 

0.13 

0.14 

^29 

0.21 

0.09 

0.11 

0.09 

0.13 

0.11 

NORSAR  RMS  Lg" 

0.27 

IQgll 

0.08 

(assuming  yields  are  subject  to  10%  S.E.  only) 

m, 

0.24 

0.16 

0.13 

0.15 

0.16 

CM 

_ 1 

0.21 

0.13 

0.11 

0.13 

0.14 

^2.9 

0.18 

0.10 

0.08 

0.10 

0.11 

NORSAR  RMS  Lg 

0.27 

0.13 

0.08 

0.09 

0.11 

*)  19  Semipalatinsk  shots  in  Table  3  used  as  calibration  events. 

**)  9  Semipalatinsk  events  with  RMS  reported  in  Ringdal  (1990)  and  Ringdal  and  Marshall 
(1989). 
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mbrYield  Relationship  at  KTS 


Yield(KT)  (Bocharov  et  al.,  1989) 

OWLS  (uncertain  X  &  Y):  5=0.80(0.024),  1=4.29(0.039),  19.  data  used, 

95%  error  in  mb  at  1,1 0,50,1 00,1 50KT;  0.24, 0.15,  0.13,  0.15,  0.16, 

95%  factor  in  yield  at  1,1 0,50,1 00,1 50KT;  3.96,  2.44,  2.13,  2.38,  2.57 

OWLS  (precise  X  assumed):  5=0.82(0.032),  1=4.26(0.052) 

Standard  LS:  5=0.81(0.031),  1=4.28(0.050) 

10%  s.e.  in  yields  assumed;  simpie  network-averaged  mb  used 


Scatter  Plot  of  Inferred  Parameters 


95%  contidcnco  Intcrv.nl  o(  slope:  0.00+/-0.051 
95%  confidence  interval  of  Intercept:  4.29+/-0.002 
[97.5%  quantile  of  t(17.  O.o.F.),  2.110,  used) 

Figure  11.  Regressing  the  simple  network-averaged  nit,  {i.e  .  m, )  on  me  i9  Soviei-pubiished 
yields.  The  yields  are  assumed  to  be  subject  to  10%  standard  errors  The  un,-;enain!ies  in  the  mj,s 
and  the  yields  are  taken  into  account  through  800  bootstrap  re.samplings  The  darkened  bundle  is 
actually  the  collection  ol  all  800  regressions,  each  produced  by  a  possible  reali/  iiion  ol  19  peiiurbed 
{nil,,  yield)  p.iirs,  I  he  95%  coniKleiic.e  band  (shown  as  2  curves  attriind  the  d.irl  .  iii'd  bundle)  i .  ii.i/ 
rower  neat  the  ceiilioid  anrl  wider  Inwards  Irolh  ends,  as  expected  I  he  individual  '.I'r/,.  conlidein^e 
intervals  ol  the  two  interred  parameters  {i.e.,  the  slope  and  the  intercept  ot  the  calibration  curve)  are 
shown  with  the  dashed  line  in  the  scatter  plot  (bottom).  Note  that  the  dashed  lectangle  is  not  the 
joint  90%  confidence  interval,  however,  due  to  the  highly  correlated  nature  ol  itre  iwo  parameters 
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800-bootstrap  Intercepts  mb(WWSSN,TG,Pmax) 


mb:Yield  Relationship  at  KTS 


Yield{KT)  (Bocharov  et  al.,  1989) 

OWLS  (uncertain  X  &  Y):  5=0.80(0.021),  1=4.32(0.033),  19.  data  used, 

95%  error  In  mb  at  1,10,50,1 00,1 50KT:  0.21, 0.13,  O.ll,  0.13,  0.14, 

95%  factor  In  yield  at  1,1 0,50,1 00,1 50KT:  3.45,  2.09, 1.86,  2.07,  2.23 

OWLS  (precise  X  assumed):  5=0.82(0.032),  1=4.28(0.053) 

Standard  LS:  5=0.80(0.031),  1=4.31(0.049) 

10%  s.o.  In  yields  assumed;  GLM  station  terms  removed 


Scatter  Plot  of  Inferred  Parameters 


95%  conlldence  Interval  o(  slope:  0.80+/m.044 
95%  confidence  Interval  of  intercept:  4.32+/-0.070 
(97.5%  quantile  of  f(l7.  D.o.F.),  2.110,  used) 


Figure  12.  Snnio  .i'-,  f  igure  11  except  Iho  f»i,s  arc  tliosc  wilh  slalion  correclions  .ippheif,  namely 
llie  nh  ;  Note  that  llie  yields  are  assumed  to  have  10%  standard  error  as  in  I  igure  11,  and  the 
reduction  in  the  scalier  is  due  to  the  bettor  precision  of 
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800-bootstrap  Intercepts  mb(WWSSN,TG,Pmax) 


mb:Yiolcl  Relationship  at  KTS 


OWLS  (uncertain  X  &  Y);  3=0.81(0.018),  1=4.28(0.027),  19.  data  used, 

95%  error  In  mb  at  1,1 0,50,1 00,1 50KT:  0.18, 0.10, 0.08, 0.10,  0.11, 

95%  factor  in  yield  at  1,1 0,50,1 00,1 50KT:  2.83, 1.74, 1.60, 1.77, 1.91 

OWLS  (precise  X  assumed):  3=0.80(0.030),  1=4.29(0.047) 

Standard  LS:  3=0.81(0.029),  1=4.27(0.046) 

10%  s.e.  in  yields  assumed;  GLM  station  &  near-source  effects  removed 


Scatter  Plot  of  Inferred  Parameters 
4.50-1 - - - 


4.40  H 


800-bootstrap  Slopes 


95%  conlidoiicc  interval  ol  slope;  0.01+/-0.037 
95%  conlidcnco  Interval  of  Intercept:  4.20+/-0.057 
[97.5%  quantile  o(  t(17.  D.o.F.),  2.110,  used] 


Figure  13.  Same  as  Figure  11  except  the  m^s  are  those  with  both  (he  slatio.-!  correciions  and 
near-source  focusing  terms  applied,  namely  (he  The  improved  precision  has  direct  impact 
on  the  regression,  as  compared  to  f  igures  11  and  12.  The  yield  of  future  Scmipatatinsk  e>plosions 
can  be  reliably  predicted  using  this  /)l,,9;yield  calibration  curve  for  a  precision  smuiar  to  .vtial  RMS 
I ,,  could  provide  {cf.  Figure  H) 
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mb:Yield  Relationship  at  KTS 


OWLS  (uncertain  X  &  Y);  3=0.71(0.028),  1=4.54(0.050),  9.  data  used, 

95%  error  In  mb  at  1,1 0,50,1 00,1 50KT:  0.27, 0.13,  0.08,  0.09,  0.10, 

95%  factor  in  yield  at  1, 10, 50,1 00,1 50KT:  5.74,  2.39, 1 .65, 1 .80, 1.98 

OWLS  (precise  X  assumed):  3=0.71(0.035),  1=4.53(0.063) 

Standard  LS:  3=0.71(0.033),  (=4.53(0.058) 

10%  s.c.  in  yields  assumed 

Scatter  Plot  of  Inferred  Parameters 


9fj7o  confidence  interval  of  slope:  0.71+/-0.067 
957o  conlidonce  interval  of  intercept:  4.54+/-0.110 
(97.5%  quanfile  of  1(7.  D.o.F.),  2.365,  used) 


Figure  14.  Regressing  9  RMS  Lg  values  (Ringdal,  1990)  on  Soviet-published  yields  RMS  Lg 
recorded  at  NORSAR  has  very  high  S/N  ratio  and  hence  very  stable  source  measure  lor  Semipala- 
tinsk  explosions  above  10KT.  There  is  no  calibration  data  below  lOKT  and  hence  the  extrapolation 
lor  future  oveni.s  with  RMS  Lg  would  have  inherently  larger  uncertainly,  as  illustrated  by  the  much 
wider  95%  conlidence  band  Thus  either  teleseismic  records  based  on  P  phases  or  Lg  measured  at 
in  state  stations  must  be  used  lor  low-yiekj  events. 


l. 6  GEOPHYSICAL  INTERPRETATION  OF  OUR  NEAR-SOURCE  CORRECTIONS 

If  we  remove  the  globally-averaged  source  sizes  from  the  station-corrected  mag¬ 
nitudes,  all  three  test  sites  would  exhibit  different  azimuthal  and  radial  amplitude  vari¬ 
ations  (Figure  5):  Degelen  and  Murzhik  events  are  systematically  enhanced  in  the 
western  U.S.  and  reduced  in  eastern  U.S.,  whereas  Balapan  events  are  all  reduced  in 
the  whole  U.S.  Degelen  events  are  reduced  in  Indonesia  and  southeast  Africa, 
whereas  Balapan  events  are  enhanced  in  these  regions.  Murzhik  events  are  reduced 
in  Scandinavia,  but  Balapan  and  Degelen  events  get  enhanced  there.  Such  hignly 
direction-dependent,  distance-dependent,  and  site-dependent  patterns  of  the  ampli¬ 
tude  fluctuation  could  be  a  diagnostic  for  the  path  effects  in  the  proximity  of  the  test 
sites.  Back  projections  (e.g. ,  Lynnes  and  Lay,  1990)  of  the  m^,  residuals  onto  the 
upper  mantle  and  the  lower  crust  reveal  that  similar  residuals  come  into  align¬ 
ment  in  several  regions  partitioned  by  known  geological  features  (Figure  15).  Murzhik 
events  recorded  in  the  western  U.S.  and  in  northeast  Asia,  Degelen  events  in  the 
western  U.S.,  and  SW  Balapan  events  at  western  European  stations  must  pass 
through  the  area  between  Chinrau  fault  and  Chingiz-Kalba  shear  zone.  All  these 
paths  show  positive  residuals.  The  north  of  Chinrau  fault  might  have  smaller  P„ 
velocity  and  higher  heat  flow  (Bonham  etal.,  1980;  Leith,  1987a,  1987b)  and  has 
negative  mean  m^,  residuals  on  the  back  projections.  Paths  from  NE  Balapan  to  North 
America  and  many  continental  European  stations  must  cross  this  area  or  even  travel 
along  the  Chinrau  fault  before  entering  deeper  mantle,  and  hence  the  complexity  in 
the  waveforms  is  inevitable.  It  seems  that  the  mean  mi^-Lg  separation  of  0.14±0.02 

m. u.  {e.g.,  Ringdal  and  Hokland,  1987;  Ringdal  and  Marshall,  1989;  Richards  etal., 
1990,  Jih  and  Wagner,  1990)  between  the  NE  and  SW  subregions  of  Balapan  could 
be  due  in  part  to  the  path  effects  —  in  addition  to  the  difference  of  source  medium 
postulated  previously  by  Marshall  et  al.  (1984).  Path  effects  can  also  explain  why  the 
SW  Balapan  waveforms  tend  to  be  more  complex  at  YKA  than  those  recorded  at 
WRA,  EKA,  and  GBA  arrays  (Jih  and  Wagner,  1991). 
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The  initial  P  waves  from  the  three  adjacent  test  sites  have  virtually  the  same 
incident  angle  at  each  teleseismic  station,  and  anything  in  common  across  all  events 
(such  as  the  crustal  amplification  as  well  as  the  upper  mantle  attenuation  underneath 
the  receiver)  would  have  been  lumped  into  the  constant  station  term.  Thus  the  station 
residuals  averaged  over  all  events  from  the  same  test  site  would  correlate  very  little 
with  the  receiver.  Instead,  they  should  reveal  more  site-dependent  information  about 
the  focusing/defocusing  pattern  underneath  E.  Kazakhstan  (Figure  15). 

The  largest  and  prominent  fault  in  the  region  is  the  southeast-trending  Chingiz 
right-lateral  strike-slip  fault  that  passes  about  10  km  southwest  of  Degelen  Mountain 
and  right  across  the  Murzhik  test  area  (Rodean,  1979;  Bonham  etal.,  1980;  Leith, 
1987b).  Soviets  reported  that  this  fault  has  a  very  steep  dip,  which  is  consistent  with 
its  linear  expression  over  large  distance  as  seen  on  Landsat  imagery  (Bonham  et  al. , 
1980).  A  distinct  fault-line  scarp  is  developed  along  much  of  the  oldest  metamorphic 
rocks.  Chingiz  Fault  extends  for  a  total  length  of  about  700  km.  Soviet  reports  postu¬ 
late  that  this  fault  extends  down  to  the  boundary  of  the  granite  layer  of  the  crust  and 
possibly  into  the  upper  mantle.  For  Murzhik  explosions,  the  propagation  of  P„  and  Lg 
waves  could  be  affected  by  this  fault  significantly,  which  results  in  a  radiation  pattern 
such  as  we  are  observing.  More  specifically,  the  rays  towards  NW  direction  could  be 
reflected  or  diffracted  to  other  quadrants,  due  to  its  post-critical  incidence  angles. 
Such  relatively  distant  crustal  structure  should  have  little  impact  on  the  first  P  waves 
of  Balapan  explosions  at  teleseismic  distances,  however. 
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GRIDDED  LS-AVERAGED  mb[Pmax]  RESIDUALS  OF  E.  KAZAKH  SHOTS 
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Figure  15.  Back-projected  residuals  averaged  over  a  grid  of  every  0  2  d-jqree  by  0  2  degree 
underneath  Soviet’s  Semipalatinsk  nuclear  test  sites.  The  averaged  residual  pattern  show  sonie 
association  with  local  geological  teatures  if  the  residuals  are  projccied  to  depths  down  to  the  lower 
crust  and  the  upper  mantle.  Thus  even  the  teleseismic  recordings  might  revi.Ml  some  uselul  infor 
rnalion  about  the  path  ettects  on  the  propagation  of  regional  phases  such  as  l'„  amt  I  ,. 
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1.7  DISCUSSION  AND  CONCLUSIONS 


The  new  magnitude  determination  scheme  (Equation  [3])  presented  in  this  study 
significantly  reduces  the  fluctuational  variation  across  the  recording  stations.  It  is 
shown  that,  by  applying  this  scheme  to  worldwide  explosions,  it  is  possible  to  have  a 
consistent  base  line  in  estimating  the  absolute  magnitudes  (which  is  crucial  in  estimat¬ 
ing  the  test  site  bias)  while  the  precision  in  the  resulting  network  m^j  values  can  be 
maintained  as  well  as  could  be  achieved  by  the  single-test-site  approach.  The  stan¬ 
dard  error  in  most  rtT2  9  values  is  0.02  m.u.,  about  the  same  as  that  for  RMS  Lg 
inferred  from  in-country  regional  network  recordings  reported  by  Israelson  (1991)  and 
Hansen  et  al.  (1990).  Most  Murzhik  events  show  nearly  identical  residual  patterns, 
suggesting  a  common  focusing/defocusing  structure.  Further  partitioning  of  Balapan 
test  site  seems  necessary,  however. 

The  most  detailed  description  of  the  wave  propagation  model  would  naturally  sug¬ 
gest  that  yet  another  term  could  be  added  into  Equation  [3]  to  count  for  the  source- 
region  attenuation.  So  far  such  source  region  bias  in  has  always  been  inferred 
with  other  information  such  as  the  yields  (as  in  this  study).  Mg  (e.g. ,  Evernden  and 
Marsh,  1987)  or  P„  velocity  (e.g.,  Marshall  etal.,  1979)  etc..  The  hypothesis  that 
such  attenuation  differential  could  be  directly  discerned  from  alone  by  further 
improving  Equation  [3]  is  worth  testing. 

Digital  signals  recorded  on  seismic  instruments  at  regional  distance  will  be  critical 
for  monitoring  low  yield  explosions  below  10  kt.  Obviously  the  regression  routines 
developed  in  this  study  can  well  be  applied  to  other  source  measures  which  are  based 
on  regional  phases.  On  the  other  hand,  teleseismic  data  such  as  the  WWSSN  data 
used  in  this  study  still  carry  invaluable  information  that  is  worthy  of  further  exploitation, 
beyond  simply  calculating  the  “unified  yield". 

Throughout  this  study,  our  emphasis  has  been  to  reveal  the  site-dependent 
characteristics  from  the  observations  exclusively  so  that  in  the  future  the  inferred 
results  can  be  critically  examined  and  compared  with  those  derived  by  other  means. 


We  suggest  that  the  follow-up  research  be  accompanied  by  well-constrained  forward 
modeling  studies  using  realistic  structures.  The  upgraded  linear  finite-difference  code 
which  incorporates  boundary  conditions  for  topographical  tree-surface  of  arbitrary 
shape  (Jih  et  al. ,  1988)  in  addition  to  the  “strain  filter”  and  “marching  grid"  features  as 
outlined  in  Jih  et  al.  (1989)  can  be  utilized  in  the  future  to  improve  our  understanding 
of  the  fundamental  issues  of  seismic  energy  partitioning  on  the  focal  sphere  as  well  as 
their  implications  for  yield  determination. 
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APPENDIX  A 


GEOTECH’S  MAXIMUM-LIKELIHOOD  NETWORK  GLM91A 


Short-period  WWSSN  vertical  recordings  (SPZ)  of  body  waves  from  Soviet 
nuclear  explosions  detonated  at  the  Semipalatinsk  Test  Site,  Eastern  Kazakhstan, 
USSR,  are  being  measured  and  added  to  our  database  to  determine  the  optimal  net¬ 
work  magnitudes  using  the  maximum-likelihood  estimator  (MLE),  which  accounts  for 
the  effects  of  data  censoring  due  to  clipping  and  to  noise  (Blandford  and  Shumway, 
1982;  Jih  and  Shumway,  1989).  As  of  now,  our  WWSSN  database  has  been 
expanded  to  192  events  (totaling  515  usable  “a”,  “b”,  and  "max”  event  phases)  from 
a  variety  of  test  sites.  Only  the  stations  at  teleseismic  distance  (20  to  95  degrees) 
which  recorded  7  or  more  good  signals  were  used  in  the  network  determination. 
Although  we  have  also  included  some  measurements  made  off  LRSM  tapes  and 
CDSN  recordings,  most  of  those  data  failed  to  meet  the  criteria  aforementioned.  The 
12170  good  signals,  8047  noise  measurements,  and  1330  clipped  recordings  yield  a 
^MLE  of  0.300. 

The  192  events  in  Table  A.2  are  grouped  by  test  sites.  Ill  events  were  meas¬ 
ured  before  1/1/90  under  various  contracts  during  the  past  decade.  25  Balapan  events 
(with  prefix  “SAP”)  were  based  on  the  raw  WWSSN  station  magnitudes  distributed  by 
DARPA  in  1988.  The  three  numbers  under  the  column  “#  of  signals”  represent  the 
number  of  signals,  noise,  and  clips  associated  with  the  P^ax  phase  of  each  event. 
S.E.M.  is  the  “standard  error  in  the  mean"  of  the  event  magnitudes.  Except  for  the 
U.S.  and  French  Sahara  explosions  which  have  specific  code  names,  all  the  remaining 
events  are  identified  with  the  dates  and  abbreviated  test  site  codes  shown  below: 


Table  A.1. 

Geotech’s  Database 

Code 

Number  of  Events 

Nuclear  Test  Site 

1/1/90 

7/15/91 

19 

37 

Nevada  Test  Site,  U.S.A. 

6 

6 

Outside  Nevada  Test  Site,  U.S.A. 

3 

3 

Amchitka  Island,  Aleutians,  U.S.A. 

A2G 

11 

11 

/^gir,  U.S.S.R. 

PNE 

1 

2 

“PNE",  U.S.S.R. 

MEK 

0 

14 

Murzhik  (Konystan),  E.  Kazakh,  U.S.S.R. 

DEK 

9 

21 

Degelen  Mountain,  E.  Kazakh,  U.S.S.R. 

SEK 

12 

22 

Balapan  (Shagan  River),  E.  Kazakh,  U.S.S.R. 

SAP 

0 

25 

Balapan  (Shagan  River),  E.  Kazakh,  U.S.S.R. 

NNZ 

18 

18 

Northern  Novaya  Zemlya,  U.S.S.R. 

SNZ 

6 

6 

Southern  Novaya  Zemlya,  U.S.S.R. 

- 

9 

9 

Ahaggar,  French  Sahara 

TU 

11 

11 

Tuamoto  Islands,  France 

RAJ 

1 

1 

Rajasthan,  India 

CH 

6 

6 

Lop  Nor,  Sinkiang,  China 
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Table  A.2.  Geotech’s  Maximum-Likelihood  Network 


Event 


ALMENDRO 


BANEBERRY 


BENHAM 


BILBY 


BOURBON 

BOXCAR 


CALABASH 


CAMBRIC 


CARPETBAG 

CHANCELLOR 

CHARTREUSE 


CHATEAUGAY 

COMMODORE 

CORDUROY 

DISCUSTHROWER 

DURYEA 

FLASK 

GREELEY 

HALFBEAK 

HANDCAR 


#  of  Signals 
26  0  2 


14  30  0 


42  1  7 


36  3  0 


18  31  0 


14  35  0 


31  16  1 


17  28  2 
31  5  1 

18  14  0 
12  39  1 
23  29  0 
36  8  0 
49  2  2 
43  2  2 
16  33  0 


S.E.M. 

.  .  j 

i^max ) 

mb(Pb) 

mb(Pa) 

0.057 

6.229 

6.021 

5.732 

0.045 

4.869 

4.547 

4.373 

0.042 

6.392 

6.140 

5.829 

0.048 

5.706 

5.453 

5.201 

0.043 

4.931 

4.751 

4.621 

0.050 

6.443 

6.220 

5.887 

0.041 

5.551 

5.357 

5.180 

0.043 

4.576 

4.310 

4.012 

0.045 

5.806 

5.585 

5.352 

0.058 

5.360 

5.201 

4.924 

0.043 

5.260 

5.029 

4.909 

0,044 

5.080 

4.898 

4.509 

0.049 

5.794 

5.585 

5.361 

0,053 

5.324 

5.131 

5.013 

0.042 

4.677 

4.451 

0.042 

5.049 

4.874 

4.723 

0.045 

5.509 

5.221 

5.038 

0.041 

6.340 

6.143 

5.909 

0.044 

6.113 

5.811 

5.583 

0.043 

4.650 

4.516 

4.345 
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Table  A.2.  Geotech’s  Maximum-Likelihood  Network  (Continued) 


Event 


HANDLEY 


HARZER 


KANKAKEE 


KNICKERBOCKER 


MAST 


MINIATA 


NASH 


PALANQUIN 


PILEDRIVER 


PURSE 


REX 


SCAUP 


SCHOONER 


SCOTCH 


SCROLL 


STARWORT 


STILTON 


CANNIKIN 


LONGSHOT 


MILROW 


#  of  Signals 

S.E.M. 

41  1  1 

0.046 

6.519 

31  5  1  ' 

0.049 

5.549 

24  27  0 

0.042 

4.875 

28  21  0 

0.043 

5.253 

29  1  0 

0.055 

6.040 

37  7  0 

0.045 

5.491 

31  21  0 

0.042 

5.166 

2  0  0 

0.212 

3.942 

40  12  2 

0.041 

5.480 

9  0  0 

0.100 

5.880 

16  35  1 

0.042 

4.778 

2  1  0 

0.173 

4.625 

7  9  0 

0.075 

4.389 

38  8  1 

0.044 

5.643 

2  0  0 

0.212 

4.077 

21  6  0 

0.058 

5.474 

7  0  0 

0.114 

5.839 

49  0  20 

0.036 

6.957 

71  4  3 

0.034 

5.873 

52  0  4 

0.040 

6.544 

6.345 


5.327 


4.624 


4.976 


5.800 


5.176 


4.939 


4.305 


4.371 


5.386 


3.642 


5.162 


5.663 


6.710 


5.494 


6.245 


6.100 


5.031 


4.390 


4.778 


5.465 


4.908 


4.789 


4.996 


5.296 


3.952 


4.247 


3.869 


5.133 


4.937 


5.455 


6.463 


5.137 


6.000 


Table  A.2.  Geotech’s  Maximum-Likelihood  Network  (Continued) 


Event 

#  of  Signals 

S.E.M. 

{Pmax ) 

FAULTLESS 

47  1  3 

0.042 

6.497 

GASBUGGY 

11  37  0 

0.043 

4.690 

RIOBLANCO 

15  20  0 

0.051 

4.831 

RULISON 

9  37  0 

0.044 

4.595 

SALMON 

6  33  0 

0.048 

4.200 

SHOAL 

16  27  0 

0.046 

4.776 

AZG22APR66 

3  10  0 

0.083 

4.225 

AZG01JUL68 

44  10  3 

0.040 

5.542 

AZG22DEC71 

12  0  2 

0.080 

6.181 

AZG25APR75 

1  16  0 

0.073 

3.986 

AZG29JUL76 

41  5  7 

0.041 

5.877 

AZG30SEP77 

21  30  1 

0.042 

4.855 

AZG170CT78 

7  0  5 

0.087 

6.108 

AZG18DEC78 

9  0  3 

0.087 

6.155 

AZG17JAN79 

10  0  4 

0.080 

6.170 

AZG14JUL79 

10  0  1 

0.091 

5.725 

AZG240CT79 

3  0  6 

0.100 

5.942 

mbiPa) 


6.193 


4.438 


4.568 


4.287 


3.989 


4.497 


5.869 


4.197 


4.127 


4.161 


3.484 


4.346 


4.144 


5.245 


5.845 


3.919 


4.932 


5.490 


5.133 


4.092 


5.294 


5.406 


5.524 

4.866 

4.865 


Table  A.2.  Geotech's  Maximum-Likelihood  Network  (Continued) 


Event 


PNE21MAY68 


PNE29AUG74 


KON18DEC66 


KON16SEP67 


KON22SEP67 


KON22NOV67 


KON31MAY69 


KON28DEC69 


KON21JUL70 


KON04NOV70 


KON06JUN71 


KON19JUN71 


KON09OCT71 


K0N210CT71 


KON26AUG72 


KON02SEP72 


DEK21NOV65 


#  of  Signals 


41  9  1 


S.E.M. 


^b(^max ) 


mbiPb)  IT^biPa) 


0.042 


0.045 


5.301 


4.753 


5.089 


4.433 


4.858 


4.041 


55  8  1 

0.038 

5.726 

36  29  2 

0.037 

5.086 

35  31  1 

0.037 

5.013 

7  64  0 

0.036 

4.317 

30  31  0 

0.038 

4.990 

45  9  3 

0.040 

5.652 

38  21  1 

0.039 

5.178 

38  22  1 

0.038 

5.242 

38  12  2 

0.042 

5.321 

41  13  0 

0.041 

5.297 

27  12  3 

0.046 

5.165 

32  9  0 

0.047 

5.359 

29  15  2 

0.044 

5.155 

15  29  0 

0.045 

4.602 

4.757 


4.001 


4.775 


5.468 


4.933 


5.053 


5.119 


5.076 


4.977 


5.139 


4.934 


4.330 


5.280 

4.558 


4.447 


4.398 


5.192 


4.592 


4.844 


4.793 


4.783 


4.742 


4.795 


.606 


^9 


Table  A.2.  Geotech's  Maximum-Likelihood  Network  (Continued) 


Event 


#  of  Signals 


DEK190CT66 


DEK26FEB67 


DEK29SEP68 


DEK23JUL69 


DEK11SEP69 


DEK22MAR71 


DEK25APR71 


DEK30DEC71 


DEK23MAR72 


DEK16AUG72 


DEK10DEC72 


DEK29MAR77 


DEK30JUL77 


DEK26MAR78 

DEK22APR78 


DEK28JUL78 


DEK22MAY80 


51  10  5 


48  9  6 


50  8  6 


38  21  1 


19  39  0 


43  14  3 


37  5  0 


163  0 


28  17  0 


23  23  1 


30  7  5 


25  14  0 


21  16  0 
25  6  0 
21  9  0 
36  9  6 
36  23  1 


S.E.M. 

^bi^max ) 

0.037 

5.539 

0.038 

5.833 

0.038 

5.631 

0.039 

5.183 

0.039 

4.603 

0.039 

5.519 

0.046 

5.783 

0.069 

5.553 

0.045 

4.979 

0.044 

4.905 

0.046 

5.542 

0.048 

5.004 

0.049 

4.877 

0.054 

5.507 

0.055 

5.020 

0.042 

5.524 

0.039 

5.129 

mbiPt)  nibiPa) 


5.439 


4.943 


4.265 


5.347 


5.594 


5.377 


4.747 


4.650 


5.340 


4.723 


4.630 


5.284 

4.771 

5.313 


4.926 


5.112 


5.368 


5.138 


4.628 


4.013 


5.052 


5.331 


5.020 


4.380 


4.361 


4.990 


4.329 


4.230 


4.963 


4.480 

5.002 


4.671 


Table  A.2.  Geotech’s  Maximum-Likelihood  Network  (Continued) 


Event 

#  Of  Signals 

S.E.M. 

(^max ) 

^b(Pb) 

SEK30JUN71 

31  19  1 

0.042 

5.062 

4.794 

SEK10FEB72 

34  8  2 

0.045 

5.319 

5.089 

SEK02NOV72 

42  1  15 

0.039 

6.185 

5.944 

SEK10DEC72 

44  2  11 

0.040 

6.020 

5.794 

SEK23JUL73 

53  1  1 

0.041 

6.191 

6.006 

SEK14DEC73 

49  8  6 

0.038 

5.760 

5.564 

SEK27APR75 

18  1  1 

0.067 

5.494 

5.254 

SEK04JUL76 

38  0  5 

0.046 

5.848 

5.601 

SEK07DEC76 

17  2  1 

0.067 

5.615 

5.420 

SEK11JUN78 

170  1 

0.071 

5.879 

5.572 

SEK15SEP78 

37  1  6 

0.045 

5.851 

5.698 

SEK23JUN79 

40  3  3 

0.044 

6.060 

5.860 

SEK04AUG79 

40  5  20 

0.037 

6.093 

5.861 

SEK280CT79 

44  5  13 

0.038 

5.946 

5.706 

SEK23DEC79 

41  3  17 

0.038 

6.145 

5.894 

SEK14SEP80 

34  5  6 

0.045 

6.033 

5.771 

SEK180CT81 

41  4  7 

0.042 

5.979 

5.754 

SEK26MAY84 

30  0  3 

0.052 

6.002 

5.915 

j  SEK14SEP88 

25  0  1 

0.059 

6.034 

5.762 

mb(Pa) 
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Table  A.2.  Geotech's  Maximum-Likelihood  Network  (Continued) 


Event 


SAF23NOV76 


SAF29AUG78 


SAF29NOV78 


SAF07JUL79 


SAF18AUG79 


SAF02DEC79 


SAF12OCT80 


SAF14DEC80 


SAF27DEC80 


SAF22APR81 


SAF13SEP81 


SAF27DEC81 


SAF25APR82 


SAF04JUL82 


SAF05DEC82 


SAF12JUN83 


SAF06OCT83 


SAF260CT83 

SAF25APR84 


22  0  0 


16  0  0 


28  0  0 


30  0  0 


28  0  0 


15  0  0 


23  0  0 


29  0  0 


24  0  0 


25  0  0 


1700 


23  0  0 


14  0  0 


21  0  0 


26  0  0 


16  0  0 
25  0  0 


18  0  0 
21  0  0 


S.E.M. 

(^max ) 

0.064 

5.626 

0.075 

5.905 

0.057 

5.880 

0.055 

5.799 

0.057 

6.087 

0.078 

5.874 

0.063 

5.828 

0.056 

5.911 

0.061 

5.896 

0.060 

5.865 

0.073 

6.024 

0.063 

6.207 

0.080 

5.944 

0.066 

6.089 

0.059 

6.093 

0.075 

5.943 

0.060 

5.942 

0.071 

0.066 
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Table  A.2.  Geotech's  Maximum-Likelihood  Network  (Continued) 


Event 

#  of  Signals 

S.E.M. 

(^max ) 

mbiPb) 

mbiPa) 

SAF270CT84 

19  0  0 

0.069 

6.150 

SAF02DEC84  22  0  0  0.064 

SAF16DEC84  15  0  0  0.078 


SAF28DEC84 


19  0  0 


0.069 


5.916 


SAF15JUN85 

15  0  0 

0.078 

6.069 

NNZ270CT66 

56  0  14 

0.036 

6.447 

6.305 

6.075 

NNZ210CT67 

53  5  3 

0.038 

5.783 

5.611 

5.424 

NNZ07NOV68 

59  1  5 

0.037 

6.042 

5.847 

5.602 

NNZ140CT69  59  2  7  0.036  6.144  5.972  5.778 

NNZ14OCT70  35  0  22  0.040  6.820  6.640  6.436 

NNZ27SbP71  23  0  21  0.045  6.629  6.487  6.276 

NNZ28AUG72  32  0  1 1  0.046  6.383  6.261  6.008 

NNZ12SEP73  23  0  21  0.045  6.770  6.677  6.356 

NNZ29AUG74  25  0  18  0.046  6.583  6.402  6.141 

I  NNZ210CT75  23  0  17  0.048  6.548  6.344  6.110 

NNZ23AUG75  27  0  12  0.048  6.495  6.376  6.128 

NNZ20OCT76  25  34  0  0.039  4.680  4.369  4.056 

j  NNZ01SEP77  25  2  2  0.056  5.572  5.433  5.126 

— 

NNZ10AUG78  39  3  18  0.039  5.867  5.637  5.414 

NNZ11OCT80  42  4  6  0.042  5.674  5.460  5.202 
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Table  A.2.  Geotech’s  Maximum-Ukelihood  Network  (Continued) 


Event 

#  of  Signals 

S.E.M. 

(^max ) 

^biPa) 

NNZ01OCT81 

43  4  5 

0.042 

5.666 

5.505 

5.251 

NNZ18AUG83 

30  4  5 

0.048 

5.721 

5.542 

5.339 

_ 

NNZ250CT84 

22  3  4 

0.056 

5.610 

5.439 

5.174 

SNZ27SEP73 

32  3  1 

0.050 

5.754 

5.518 

5.227 

SNZ270C73A 

14  0  24 

0.049 

7.082 

6.864 

6.645 

SNZ270C73B 

9  28  0 

0.049 

4.189 

4.037 

SNZ270C73C 

4  34  0 

0.049 

3.951 

3.928 

3.587 

SNZ02NOV74 

12  0  29 

0.047 

7.001 

6.784 

6.502 

SNZ180CT75 

21  0  21 

0.046 

6.838 

6.527 

6.245 

BERYL 

1160 

0.073 

5.017 

4.815 

4.455 

CORUNDON 

11  41  0 

0.042 

4.247 

3.951 

3.852 

EMERAUDE 

14  25  0 

0.048 

4.596 

4.269 

GRENAT 

32  31  1 

0.038 

4.787 

4.524 

4.332 

OPALE 

3  50  0 

0.041 

. 

3.950 

_ 1 

3.909 

3.827 

RUBIS 

45  5  0 

0.042 

5.434 

5.185 

4.863 

SAPHIR 

55  5  5 

0.037 

5.725 

5.479 

5.196 

TOURMALINE 

27  39  0 

0.037 

4.671 

4.463 

4.158 

H 

TURQUOISE 

11  53  0 

0.038 

4.258 

3.986 
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Table  A.2.  Geotech’s  Maximum-Likelihood  Network  nrif,  (Continued) 


Event 

#  of  Signals 

S.E.M. 

^bi^max ) 

mb(Pb) 

nibiPa) 

TU19FEB77 

16  27  0 

0.046 

4.665 

4.413 

TU19MAR77 

20  5  1 

0.059 

5.682 

5.475 

5.175 

TU24NOV77 

32  0  0 

0.053 

5.702 

5.437 

5.102 

TU30NOV78 

38  7  2 

0.044 

5.635 

5.251 

4.862 

TU25JUL79 

1800 

0.071 

5.917 

5.624 

5.158 

TU23MAR80 

27  14  3 

0.045 

5.394 

5.141 

4./13 

TU19JUL80 

38  2  2 

0.046 

5.559 

5.202 

4.939 

TU03DEC80 

31  10  0 

0.047 

5.371 

5.020 

4.739 

TU25JUL82 

22  13  0 

0.051 

5.252 

5.078 

4.717 

TU19APR83 

21  1  0 

0.064 

5.533 

5.228 

5.011 

TU25MAY83 

18  0  0 

0.071 

5.764 

5.446 

5.163 

RAJ18MAY74 

7  23  0 

0.055 

4.595 

4.341 

4.081 

CH22SEP69 

30  12  0 

0.046 

5.190 

4.801 

4.409 

CH270CT75 

12  24  0 

0.050 

4.655 

4.467 

4.223 

CH170CT76 

12  33  0 

0.045 

4.610 

4.298 

4.134 

CH06OCT83 

16  12  1 

0.056 

5.207 

4.997 

4.740 

CH03OCT84 

10  12  0 

0.064 

5.020 

4.769 

4.489 

CH19DEC84  3  10  0  0.083  4.424  4.077  4.101 
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APPENDIX  B 


ESTIMATES  OF  TEST  SITE  BIAS  WITH  mj,(GLM91A) 

Since  the  near-source  correction  procedure  presented  in  this  study  has  not  been 
applied  to  regions  other  than  Semipalatinsk,  some  experiments  using  the  more  com¬ 
plete  m^(GLM91A)  values  may  be  interesting.  The  event  values  in  Table  A.2  are 
corrected  for  the  station  terms,  and  hence  are  very  similar  to  the  1772,2  discussed  in 
Section  1.3.  The  only  difference  is  that  7772,2  'S  computed  as  the  (maximum-likelihood) 
averaged  012,2  across  those  and  only  those  stations  which  reported  the  amplitude, 
whereas  mb(GLM91A)  is  inverted  jointly  with  all  515  events  and  122  stations  simul¬ 
taneously,  and  hence  the  missing  stations  are  also  included  in  the  (maximum- 
likelihood)  averaging.  The  discrepancy  in  the  two  mt,  values  is  insignificant,  however 
(cf.  Tables  2  and  A.2). 

Figure  16  shows  the  results  of  regressing  m^(GLM91A)  on  the  published  yields  of 
Semipalatinsk  and  NTS  high-coupling  explosions  with  10%  S.E.  in  yields  assumed. 
U.S.  and  U.S.S.R.  have  released  the  yields  of  roughly  equally  many  events  (Springer 
and  Kinaman,  1971,  1975;  Bocharov  etal.,  1989;  Vergino,  1989).  Note  that  the  NTS 
calibration  curve  based  on  phase  of  /nb(GLM91A)  (top  of  Figure  16)  has  the 
same  slope  and  intercept  as  those  based  on  m2  9  (Figure  13).  Also  note  that  the  NTS 
curve  based  on  the  first  arrivals  {i.e. ,  P^  phase)  has  a  smaller  slope,  and  there  seems 
to  be  a  lot  of  scatter  for  the  low  yields.  Blandford  (written  communication,  1991) 
pointed  out  that  this  might  occur  if  the  low-yield  mij{Pg)  were  biased  high  and  had  a 
lot  of  scatter  due  to  the  noise  wavelets  interfering  with  signal  wavelets. 


-62- 


mb.-YicId  Hclationship;  KTS  vs.  NTS 


Figure  16.  mb(GLM)  versus  published  yields  of  Semipalatinsk  and  NTS  high-coupling  explosions 
(with  10%  S  E.  in  yields  assumed).  U.S.  and  U  S.S.R.  have  released  the  yields  of  roughly  equally 
many  events.  Tlie  moan  KTS  NTS  test  site  bias  is  larger  than  0.35  m.u.  if  or  is  used,  and 
the  larger  scatter  in  NTS  calibration  curve  is  ignored.  The  bias  would  be  about  0  15-0  2  m  u  if  only 
granitic  shots  at  NTS  are  used. 


Table  B.1  lists  the  mean  KTS-NTS  test  site  bias  at  three  different  yield  levels. 
Results  based  on  "b”  and  “max”  phases  suggest  that  the  bias  is  yield-dependent  and 
appears  to  be  larger  at  the  lower  yield  end.  The  “mean”  KTS-NTS  m^,  bias  is  slightly 
larger  than  0.35  m.u.,  although  Figure  16  also  indicates  that  the  NTS  granite  events 
would  seem  to  be  about  0.15-0.2  m.u.  below  the  Semipalatinsk  curve. 


Table  B.1.  Estimated  Test  Site  Bias  (from  Nuttli’s  earlier  studies) 

Test  Sites 

Magnitude 

Description 

lOkt 

lOOkt 

150kt 

Balapan  -  NTS 

^bCISC).  Lg 

Nuttli  (1987) 

0.35 

0.35 

0.35 

Degelen  -  NTS 

mb  (ISC).  Lg 

Nuttli  (1987) 

0.58 

0.58 

0.58 

Estimated  Test  Site  Bias  (from  this  study)  | 

Test  Sites 

Magnitudes 

Description 

lOkt 

lOOkt 

KTS*  -  NTS" 

) 

Marshall 

1 

0.55 

0.50 

0.48 

KTS  -  NTS 

mb(Pa) 

TG 

0.38 

0.39 

0.40 

KTS  -  NTS 

mb(Pb) 

TG 

■B 

0.47 

KTS  -  NTS 

^bi^max ) 

TG 

0.42 

KTS  -  NTS**' 

^bi^max ) 

TG,  S-Cubed*’* 

0.36 

0.36 

0.36 

*)  Combining  all  UK/AWE's  Balapan.  DegBlen,  and  Murzhik  values  as  listed  in  Vergino  (1989) 
")  UK/AWE  s  NTS  m^  as  distributed  in  1987 

**■)  Murphy  (1981):  m^(S-Cubod)  =  3  92  +  0.81  log(W)  (or  NTS  high-coupling  events 


That  NTS  granite  events  might  lie  above  typical  wet-tuff  or  rhyolite  events  on  the 
-yield  calibration  curve  can  be  further  illustrated  by  some  simple  calculations  with 
the  three  events  Dougals  (1987)  analyzed.  The  announced  yields  of  events  68061 9B, 
710630B,  and  PILEDRIVER  are  <20  kt.  <20  kt,  and  56  kt,  respectively.  Our  yield  esti¬ 
mates  {cf.  the  Pmax  calibration  curves  shown  in  Figure  16).  however,  are  17  kt,  9  kt, 
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and  100  kt,  respectively,  based  on  their  corresponding  ^^(GLMOIA.  Pmax)  of  5  282, 
5.062,  and  5.480.  (Note  that  the  discrepancy  between  ^2  9  and  m^(GLM91A)  is 
insignificant.)  At  NTS,  a  high-coupling  17-kt  shot  has  an  expected  of  4.82,  which  is 
0.46  m.u.  below  that  of  68061 9B  at  Balapan.  Likewise,  a  9-kt  high-coupling  shot  at 
NTS  would  be  expected  to  have  a  around  4.58,  about  0.48  m.u.  smaller  than  that 
of  710630B.  At  KTS,  the  expected  for  100  kt  and  56  kt  would  be  5.9  and  5.69, 
respectively;  which  are  0.42  and  0.43  m.u.  larger  than  typical  NTS  shots  at  the 
corresponding  yields.  The  KTS-NTS  bias  of  0.48  (9  kt),  0.46  (17  kt),  0.43  (56  kt),  and 
0.42  (100  kt)  resemble  that  yield  dependency  as  shown  in  Table  B.1,  as  expected. 
PILEDRIVER’s  mj,(GLM91  A),  5.480,  is  about  0.22  m.u.  larger  than  that  of  a  56-kt  shot 
at  NTS,  i.e. ,  5.26.  This  result  seems  to  match  very  well  with  Ryall’s  (1985)  inference 
of  the  attenuation  differential  between  Semipalatinsk  and  NTS  using  earthquake  data 
recorded  at  seismic  stations  in  these  two  region. 

Murphy  (1981)  points  out  that  the  "statistically  significant"  m^-yield  relationship  for 
the  wet  tuff/rhyolite  explosions  at  Pahute  Mesa  and  Yucca  Flat  is 

=  3.92 -I- 0.81  log(W)  [1] 

which  happens  to  be  about  0.36  m.u.  below  our  inferred  calibration  curve  for  historical 
Semipalatinsk  explosions.  This  could  be  simply  accidental.  Nevertheless,  an  interest¬ 
ing  speculation  can  be  offered  to  explain  the  coincidence.  It  is  not  impossible  that  the 
Soviets  are  fully  aware  of  Equation  [1]  and  the  commonly  quoted  KTS-NTS  bias  of 
0.35  m.u.  (e.g. ,  OTA,  1988).  Perhaps  the  Soviets  have  purposefully  released  a  sub¬ 
set  of  their  historical  explosions  which  would  fit  a  prescribed  curve  roughly  0.35  m.u. 
above  Equation  [1].  If  this  was  indeed  the  case,  probably  the  released  19  Eastern 
Kazakhstan  events  were  not  "fudged”  othenwise  —  although  whether  they  are  truly 
representative  of  the  whole  explosion  population  would  still  remain  open  {cf.  the  dis¬ 
cussion  in  Gray  etal.,  1990).  We  could  also  argue  that,  if  the  aforementioned  specu¬ 
lation  were  valid,  then  Geotech's  measurements  of  Soviet  events  must  correlate 
very  well  with  the  magnitudes  which  the  Soviet  seismologists  have  used  in  determining 


their  own  yields. 


Using  20  NTS  tuff/rhyolite  events,  Marshall  et  al.  (1979)’s  m2  give 

m2  =  3.71  +  0.89  log(W).  [2] 

This  is  not  significantly  different  from  our  result  for  21  NTS  tuff/rhyolite  events  (Figure 
16): 

/77ft(GLM91A)  =  3.76  +  0.86  log(W)  [3] 

Combining  our  GLM/MLE-derived  WWSSN  with  RMS  Lg  values  measured  at 
NORSAR  (Ringdal  and  Marshall,  1989),  the  mt^-Lg  residuals  for  E.  Kazakh  explosions 
show  a  strong  difference  among  these  three  test  sites  (Figure  1 7).  The  SW  subregion 
of  Balapan  test  site  excites  slightly  larger  (relative  to  Lg),  whereas  all  the  remain¬ 
ing  regions  of  E.  Kazakh  test  site  have  negative  residuals.  All  studies  of  the  intrasite 
nig-Lg  bias  of  Balapan  explosions  lead  to  a  very  consistent  estimate,  namely 
0.1410.02  m.u.  (Table  B.2  and  Jih  and  Wagner,  1991). 

As  noted  in  Appendix  A,  DARPA  distributed  the  WWSSN  station  values  of  39 
large  Balapan  explosions  furnished  by  AFTAC  in  the  spring  of  1988.  Lilwall  et  al. 
(1988)  supplemented  this  data  set  with  some  ISC  recordings  in  their  analysis,  and 
they  found  that  the  event  values  based  on  Blacknest’s  Joint  Maximum-Likelihood 
(JML)  method  are  not  significantly  different  from  those  based  on  LSMF.  We  have 
incorporated  these  AFTAC's  values  into  our  database  (cf.  pages  58-59),  with  a 
compensating  correction  for  the  different  B(A)  factors.  The  LSMF  results  of  the  39 
AFT AC-measured  Balapan  events  show  a  mean  m^j-Lg  bias  of  0.14  between  SW  and 
NE  subregions  of  Balapan  test  site  (Table  B.2). 
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Table  B.2.  Mean  m^j-Lg  of  Balapan  Explosions 


Reference 

SW 

TZ 

NE 

SW-NE 

Ringdal  and  Hokland  (1987) 

0.112±0.009(?) 

— 

-0.059±0.014(?) 

0.17 

Marshall  (1987)  +  NORSAR 

0.116±0.009(26) 

0.041  ±0.01 2(10) 

-0.042±0.014(14) 

0.16 

Ringdal  and  Marshall  (1989) 

0.05±0.007(46) 

-0.02±0.009(20) 

-0.10±0.012(30) 

0.15 

TGAL  ^  NORSAR 

0.020 1 0.01 5(20) 

-0.071  lO.OI  7(8) 

0.112(0.017(8) 

0.13 

AFTAC  +  NORSAR 

-0.012±0.010(18) 

-0.113±0.020(6) 

-0.153±0.018(7) 

0.14 

Marshall  (1987)  +  Nuttli 

-0.012±0.015(17) 

-0.060±0.026(4) 

-0.118+0.014(14) 

0.11 

Table  B.3  lists  the  mean  m^j-Lg  values  at  three  test  sites  of  Eastern  Kazakhstan. 
Murzhik  events  have  smaller  relative  mi,  excitation,  as  compared  to  Balapan  and 
Oegelen  explosions.  Since  Balapan  and  Murzhik  events  essentially  followed  the  same 
depth-yield  scaling  (Jih  and  Shumway,  1991;  Jih,  1990),  the  relatively  strong  Lg  exci¬ 
tation  at  Murzhik  could  probably  be  due  to  the  smaller  size  (and  hence  shallower 
depth  of  burial)  at  Murzhik.  or  due  to  the  source  medium  (Figure  17). 


Table  B.3.  Mean  rn^-Z.^  of  Eastern  Kazakh  Explosions 

Reference 

SR 

DM 

MK 

TGAL  +  NORSAR 

-0.030±0.014(36) 

-0.047±0.034(5) 

-0.128±0.034(3) 

AFTAC  +  NORSAR 

-0.063±0.014(31) 

_(?) 

-—(?) 
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SPATIAL  PATTERN  OF  mb-Lg  RESIDUALS  OF  E.  KAZAKH  SHOTS 


77.500  79.100 


X 

0.187 

X 

0.004 

X 

0.062 
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0.047 

o 

.0.047 

O 

4>.062 

o 

.0.094 

O 

^>.187 

Figure  17.  The  spatial  pattern  of  mi,-Lg  residuals  of  Semipalatinsk  explosions  with  TG’s 
and  RMS  Lg  values  reported  at  NORSAR.  The  residual  pattern  of  Balapan  events  strongly  indicates 
significant  difference  in  the  source  medium  across  the  Chinrau  fault  separating  the  northeastern  and 
southwestern  portion  of  the  test  site,  as  reported  by  Ringal  and  Marshall  (1989)  and  Marshall  et  al. 
(1984).  The  mean  bias  between  SW  and  NE  Balapan  is  about  0.13  m  u 


mb:  Geotech’s  WWSSN  GLM/MLE  Pmax  (GLM91A,  515  events) 
(122  stations,  each  recorded  7  signals  or  more) 

NORSAR  RMS  Lg  (Ringdal,  1990;  Ringdal  and  Marshall,  1989) 
Surface  geology;  Bonham  et  al.  (1980),  Leith  (1987). 

Balapan,  SW  region,  20  events:  mb(TG)=  mb(Lg)  +  0.020(0.015) 
Balapan,  TZ  region,  8  events:  mb(TG)=  mb(Lg)  -  0.071(0.017) 
Balapan,  NE  region,  8  events:  mb(TG)=  mb(Lg)  -  0.112(0.017) 
Balapan,  36  events:  mb(TG)=  mb(Lg)  -  0.030(0.014) 

Dogolcn,  5  events:  mb(TG)r  mb(Lg)  -  0.047(0.034) 

Murzhik,  3  events:  mb(TG)=  mb(Lg)  -  0.128(0.034) 

Sedimentary  &  volcanic  rocks 
Devonian  &  Carboniferous  rock 
'!  Granitic  rocks 
iMiflIII  Limestone 
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SECTSON  li  PROJECT  OVERVIEW 

RECENT  METHODOLOGICAL  DEVELOPMENTS  IN 
MAGNITUDE  DETERMINATION  AND  YIELD  ESTIMATION 
WITH  APPLICATIONS  TO  SEMIPALATINSK  EXPLOSIONS 


Rong-Song  Jih 

Teledyne  Geotech  Alexandria  Laboratories 
314  Montgomery  Street 
Alexandria.  VA  22314-1581 


11.1  CONTRACT  NO.:  F19628-89-C-0063.  Task  1  (expired  July  1991) 

11.2  OBJECTIVES 

The  primary  objective  is  to  study  the  technical  issues  relating  to  the  seismic  esti¬ 
mation  of  the  yield  of  remote  underground  explosions.  Our  approach  is  to  improve 
both  the  numerical  and  statistical  modeling  tools  as  much  as  we  can,  and  then  apply 
the  upgraded  numerical  tools  to  the  excitation/propagation  study  of  the  teleseismic  and 
regional  phases,  and  apply  the  upgraded  statistical  tools  to  the  magnitude  determina¬ 
tion  as  well  as  the  yield  estimation  problems  with  emphasis  on  Soviet  explosions. 


11.3  RESEARCH  ACCOMPLISHED  DURING  CONTRACT  PERIOD 

11.3.1  Upgrading  of  Unbiased  Network  Estimator 

The  body  wave  magnitudes  used  in  developing  and  applying  the  magnitude-yield 
relationship  are  network  values  which  are  some  ‘‘average’’  of  the  station  m^. 
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Previously  a  major  objective  of  magnitude-calculation  research  was  to  determine  the 
network  that  is  not  biased  by  the  sample  truncation  due  to  the  limited  range  of  the 
seismometers.  Ringdal  (1976)  introduced  the  maximum-likelihood  estimator  [MLE]  to 
correct  for  the  statistical  bias  introduced  by  data  censoring  from  non-detection.  Von 
Seggern  and  Rivers  (1978)  pointed  out  the  importance  of  accounting  for  the  data  cen¬ 
soring  due  to  signal  clipping.  Blandford  and  Shumway  (1982)  derived  the  general 
linear  model  [GLM]  in  the  presence  of  data  censoring  using  the  Expectation  Maximiza¬ 
tion  [EM]  algorithm.  They  simultaneously  estimated  event  magnitudes  and  station 
corrections  in  a  maximum-likelihood  sense.  Jih  and  Shumway  (1989)  re-examined 
and  documented  the  GLM  algorithms,  and  they  discussed  the  uncertainty  assessment 
in  the  censoring  situation.  It  is  concluded  that  in  the  multi-parameter  linear  regression 
problem  with  censored  data,  the  scaling  of  o/Vdegrees  of  freedom  still  provides  an 
extremely  good  approximation  of  the  uncertainty  associated  with  each  parameter.  In 
the  case  of  non-censoring,  such  approximation  can  be  proved  to  be  "exact”.  The 
methodological  similarities  and  differences  between  the  iterative  least  squares  [ILS] 
and  the  maximum-likelihood  estimator  [MLE]  were  also  identified  in  Jih  and  Shumway 
(1989). 

A  recent  breakthrough  in  magnitude  determination  is  the  development  of  a  pro¬ 
cedure  to  account  for  the  near-source  focusing/defocusing  effects.  Jih  and  Wagner 
(1991b)  propose  to  compute  the  new  station  magnitude  mjg  for  the  i-th  event 
recorded  at  the  j-th  station  as 

m2.9(i.j)  =  logio[A(i,j)/T(i,j)]  +  B(A(i,j))  -  S(j)  -  F(k(i),j)  [1] 

where  A(i,j)  is  the  displacement  amplitude  (in  millimicrons)  and  T(i,j)  is  the  period  (in 
seconds)  of  the  P  wave.  The  B(A)  is  the  distance-correction  term.  S{j)  is  the  station 
correction,  and  F(k(i),j)  is  the  near-source  focusing  correction  for  explosions  from  the 
k(i)-th  source  region.  Some  of  the  S  terms  for  explosions  from  any  test  site  and  the  F 
terms  for  the  Semipalatinsk  area  are  listed  in  Table  1.  A  complete  list  can  be  found  in 
Jih  and  Wagner  (1991b).  This  new  magnitude  is  called  rr72  9  to  avoid  confusion  with 
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the  m3  defined  in  Marshall  et  al.  (1979)  that  corrects  for  the  source-region  attenuation 
and  station  terms  solely  based  on  published  velocity. 


Table  1.  Receiver  and  Near-source  Corrections  for  WWSSN  Stations  (partial  listing) 


Station  Term 

Near-source  Term,  F 

Station 

Code 

S 

Balapan 

Degelen 

Murzhik 

Longitude 

Latitude 

Description 

AAE 

-0.352 

-0.387 

-0.076 

-0.127 

38.766 

9.029 

Addis  Ababa,  Ethiopia 

AAM 

0.205 

0.202 

-0.073 

-0.247 

-83.656 

42.300 

Ann  Arbor,  Michigan 

AKU 

-0.048 

0.300 

0.311 

0.212 

-18.107 

65.687 

Akureyri,  Iceland 

ANP 

-0.349 

-0.167 

0.402 

0.183 

121.517 

25.183 

Anpu,  Taiwan 

AQU 

-0.133 

-0.195 

0.039 

-0.024 

13.403 

42.354 

Aquila,  central  Italy 

ATU 

0.128 

0.191 

-0.156 

0.034 

23.717 

37.972 

Athens  Univ.,  Greece 

BAG 

-0.027 

0.076 

-0.009 

-0.103 

120.580 

16.411 

Baguio  City,  Luzon  Island 

BEC 

-0.114 

0.090 

0.058 

-0.092 

-64.681 

32.379 

Bermuda-Columbia,  Atlantic 

BKS 

0.089 

-0.014 

0.101 

0.229 

-122.235 

37.877 

Byerly,  central  California 

BLA 

0.057 

-0.233 

-0.177 

-0.299 

-80.421 

37.211 

Blacksburg,  West  Virginia 

BOZ 

0.188 

-0.325 

-0.025 

-0.180 

-111.633 

45.600 

Bozeman,  Montana 

BUL 

0.003 

0.014 

-0.266 

-0.037 

28.613 

-20.143 

Bulawayo,  Rhodesia 

CHG 

-0.140 

0.240 

0.106 

0.045 

98.977 

18.790 

Chiengmai,  southeast  Asia 

-0.178 

0.114 

0.375 

0.602 

-115.083 

67.833 

Copper  Mine,  Canada 

COL 

0.181 

0.188 

0.051 

-147.793 

64.900 

College  Outpost,  Alaska 

COP 

-0.003 

0.159 

-0.276 

IB 

55.683 

Copenhagen,  Denmark 

COR 

0.155 

0.132 

0.172 

44.586 

Corvallis,  Oregon 

CTA 

0.153 

-0.072 

-0.073 

146.254 

-20.088 

Charters  Towers,  Australia 

DAG 

0.036 

-0.052 

0.086 

-18.770 

76.770 

Danmarkshavn,  Greenland 

DAV 

_ 

DUG 

-0.320 

-0.264 

-0.053 

125.575 

7.088 

Davao,  Mindanao  Island 

0.149 

0.038 

0.371 

0.352 

-112813 

40.195 

Dugway,  Utah 

EIL 

0.004 

-0.117 

-0.228 

-0.103 

34  950 

29.550 

Eilat,  Arabic  Peninsula 

ESK 

0.048 

-0.042 

0.162 

-0.327 

-3.205 

55.317 

Eskdalemuir,  Scotland 

FLO 

0.000 

-0.294 

-0.093 

-0.446 

-90.370 

38.802 

Florissant,  eastern  Missouri 
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Testing  results  indicate  that  the  procedure  described  in  [1]  has  the  following 

advantages: 

[A]  For  79  out  of  82  Semipalatinsk  events  we  have  tested,  Equation  (1]  provides 
more  stable  measurements  across  the  whole  recording  network,  as  compared 
to  the  conventional  GLM  or  LSMF  procedure  which  only  corrects  for  the  station 
terms  [cf.  Table  2).  The  reduction  in  the  standard  deviation  of  network  from 

to  r?T2  9  could  reach  a  factor  of  3  {cf.  Table  2).  Events  which  do  not  show 
improvements  in  precision  could  have  been  detonated  in  environments  with 
different  focusing  patterns. 

[B]  The  resulting  network  values  are  not  significantly  different  from  the  GLM 
results.  Thus  if  the  mean  network  mj,  values  derived  by  GLM  or  LSMF  are 
unbiased,  so  are  the  refined  results. 

[C]  The  scatter  in  /tTj  9  versus  log(yield)  is  smaller  than  that  for  other  {cf  Table 
3). 
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Table  2.  Various  Network-Averaged  mt,  of  19  Special  Events* 


Event 

#  of  signals 

Without  correction 

Station  corrected 

Near-source  corrected 

Date 

,  a 

^2.2  •  ® 

^2.9  ■  ^ 

651121D 

48  12  1 

5.38410.035  0.271 

5.39410.024  0.188 

5.38110.019  0.152 

66021 3D 

51  2  10 

6.07310.038  0.305 

6.09210.027  0.218 

6.08810.014  0.114 

660320D 

670922M 

49  6  8 

5.84810.041  0.322 

5.86510.030  0.239 

5.85310.009  0.074 

5.029+0.017  0.125 

35  20  1 

5.03310.036  0.271 

5.04810.029  0.214 

680929D 

50  4  6 

5.61010.035  0.275 

5.64210.026  0.202 

5.64110.018  0.138 

690723D 

38  17  1 

5.17210.041  0.306 

5.20110.029  0.220 

5.18610.013  0.100 

691130B 

51  0  0 

5.95010.044  0.313 

5.97310.031  0.222 

5.94510.026  0.184 

691228M 

45  2  3 

5.66610.043  0.306 

5.66510.035  0.250 

5.66010.018  0.125 

710425D 

37  3  0 

5.76410.052  0.327 

5.79310.042  0.267 

5.82610.020  0.127 

710606M 

38  6  2 

5.32310.040  0.272 

5.34110.031  0.210 

5.31910.015  0.099 

711009M 

27  9  3 

5.16510.037  0.233 

5.18710.028  0.174 

5.13610.016  0.100 

711021M 

32  6  0 

5.34810.049  0.301 

5.38310.036  0.224 

5.34110.021  0.127 

72021  OB 

34  8  2 

5.29710.042  0.278 

5.32910.029  0.195 

5.28910.018  0.121 

72081 6D 

23  20  1 

4.90810.044  0.293 

4.93110.033  0.221 

4.92110.025  0.165 

720902M 

15  25  0 

4.61510.049  0.312 

4.63510.038  0.243 

4.60210.017  0.107 

721102B 

42  0  15 

6.17310.045  0.339 

6.18310.034  0.256 

6.16010.023  0.177 

721210B 

45  1  11 

6.00610.037  0.282 

6.01 310.029  0.223 

5.98310.022  0.169 

88091 4B 

25  0  1 

6.00410.037  0.191 

6.03210.023  0.117 

6.02610.033  0.168 

*)  19  Semipalatinsk  explosions  for  which  the  yields  were  published. 

*■)  m,  =  network  average  of  raw  without  any  correction  —  equivalent  to  ISC  bulletin  nhy  =  network 
average  of  with  GLM  station  corrections  applied.  =  network  average  of  with  both  station  terms 
and  near-source  focusing  terms  removed 


Table  3.  95%  Confidence  Factor  of  Semipalatinsk  Yield  Estimate 


(assuming  yields  are  subject  to  rounding  and  10%  S.E.) 

m^  used 

1  kt 

10  kt 

50  kt 

100  kt 

150  kt 

m, 

4.74 

2.58 

2.19 

2.42 

2.61 

77122 

3.98 

2.20 

1.89 

2.12 

2.29 

^2.9 

3.23 

1.82 

1.67 

1.88 

2.04 

Ringdal's  RMS  Lg^ 

5.74 

2.39 

1.65 

1.80 

1.98 

Israelson’s  RMS  Lg^ 

2.82 

1.81 

1.72 

1.87 

1.97 

(assuming  yields  are  subject  to  10%  S.E.  only) 

3.96 

2.44 

2.13 

2.38 

2.57 

77722 

3.45 

2.09 

1.86 

2.07 

2.23 

77522 

2.83 

1.74 

1.60 

1.77 

1.91 

Ringdal’s  RMS  Lg 

5.74 

2.39 

1.65 

1.80 

1.98 

Israelson’s  RMS  Lg 

2.70 

1.84 

1.78 

1.86 

1.95 

1)  RMS  Lg  ol  9  Semipalatinsk  events  furnished  by  Ringdal  (1990)  with  NORSAR  and  GRF  data. 

2)  RMS  Lg  of  16  Semipalatinsk  events  furnished  by  Israelson  (1991b)  with  hand-digitized  Soviet  analog 
seismograms. 


11.3.2  m^j-Yleld  Regression  Routine  with  Censored  Yields:  MLE-CY 

In  general  there  are  four  types  of  yield  data  available;  [0]  the  yield,  W,  is  known 
as  yo  kt,  [1]  W  is  left  censored,  /.e. ,  the  exact  value  of  W  is  only  known  to  be  less 
than  certain  level,  [2]  W  is  right  censored,  i.e. ,  the  exact  value  of  W  is  only  known  to 
be  larger  than  certain  level,  and  [3]  W  is  only  known  to  lie  between  two  bounds.  The 
majority  of  Soviet  yields  recently  published  by  Bocharov  et  al.  (1989)  and  Vergino 
(1989)  are  of  type  3. 
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The  problem  of  estimating  the  yield  of  an  explosion  from  the  estimated  seismic 
magnitude  has  been  handled  traditionally  using  the  linear  or  piecewise  linear  model 

X  =  a  +  P  log(W)  +  V  =  a  +  pY  +  V  [2] 

where  X  is  the  measured  magnitude  (e.g.,  mi,  or  Ms),  a  and  p  are  intercept  and  slope 
estimators,  W  is  the  yield  in  kiloton  [kt],  and  v  is  an  error  term,  v  is  assumed  to  be  a 
Gaussian  random  variable  with  mean  zero  and  standard  deviation  a.  One  may  then 
collect  a  number  of  "calibration  events”,  estimating  a  and  p  by  least  squares  using  a 
number  of  known  yields  and  measured  magnitudes.  This  classical  calibration 
approach  leads  to  predicting  a  future  log-yield  Y  at  magnitude  =  X  by  inverting  Equa¬ 
tion  [2],  i.e. ,  Y  =  (X  -  <i)/p. 

The  geometrical  interpretation  of  “regressing  X  on  Y”  is  that  the  (6,  p)  thus 
estimated  will  be  the  optimal  solution  that  minimizes  the  sum  of  the  squared  X  residu¬ 
als,  2]  (  X  -  &  -  pY  )2.  Implicitly,  here  we  have  made  an  assumption  that  the  indepen¬ 
dent  variable  Y  has  nearly  perfect  accuracy  and  precision  as  compared  to  X.  Alter¬ 
nately,  one  can  estimate  X  and  q  in  the  inverse  regression  model 

Y  =  X  +  TiX-hv'  [3] 

and  then  predict  a  future  log-yield  directly  as  ?  =  X  +  f|  X.  Likewise,  here  one  is  impli¬ 
citly  assuming  that  X  has  perfect  accuracy  and  precision,  and  hence  the  optimal  esti¬ 
mate  (X,,  fj)  is  the  one  that  minimizes  the  sum  of  squared  Y  residuals, 
^  (  Y  -  ^  -  f^X  )2.  Thus  either  the  yield  or  the  magnitude  must  be  regarded  as  an 
error-free  independent  variable  in  these  two  models.  Symbolically,  these  two  conven¬ 
tional  regression  models  are  based  on  the  following  two  extreme  assumptions: 
o(X)/a(Y)  =  oo  and  o(X)/o{Y)  =  0,  respectively. 

In  reality,  both  the  and  the  yield  measurements  are  subject  to  error.  At  NTS, 
o(/T7^)  »  o(log  yield)  could  be  a  reasonable  assumption  to  justify  the  regression  of 
on  the  yields.  However,  this  may  not  be  the  case  in  general.  Note  that  [3]  can  be 
rewritten  in  a  form  similar  to  [2]:  X  =  a'  +  P'Y  +  v"  with  the  transformations  a'  =  -  Xyrj, 
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P'  =  1/Tl. 


Elegant  maximum-likelihood  theory  can  be  derived  for  “regressing  censored  Y  on 
X”  (Jih  etal.,  1990a).  Suppose  there  are  Oq,  ni,  02.  and  n3  events  for  each  type, 
respectively.  The  conditional  likelihood  function  of  the  censored  observations 
{  Vo.  '1.  ^3)  given  the  intercept  a,  slope  p,  and  a  is 


L  {  Vo.  ti.  t2.  tg  I  a,  p.  o  )  =n  P(  Yj  =  Voj  I  a.  P.  o  )  *  ri  P(  Yj  <  tij  1  a.  p.  o  )  *  [4] 

i=i  i=i 


n  P(  Yj  >  t2j  I  a,  p.  a  )  ‘  n  P(  U  <  Yj  <  tbj  |  a,  p,  o  ) 

j=i  j=i 

and  the  log-likelihood  function  is 


In  L  (  Vo.  ti.  t2,  ta  I  a,  p.  a  )  =  -  In(27ta2)  -  Z  (Voi  -  +  [5] 


ilna)(zij)+  Eln^('Z2j)+  Zlnl<l>(Zbj)-^(Zai)] 

i=i  j=i  j=i 

where  Zj  h  (a  -1-  p  t,  -  Xj)/  pa  ;  yg,  t,.  tg,  and  13  are  the  collection  of  announced  yields. 


Solving 


ainL 

da 


=  0  implies  immediately  that  the  d  must  satisfy  the  following  neces¬ 


sary  condition; 


Z(  yoi  - 


a^  = 


i-i 


,  ^  <i>(Z2i)  ^  ^  (i)(Zbi)Zbi  -  <i)(Zaj)Za, 

°  j5  O(Zbj)  -  0(Z3,) 


[6] 


Solving 


ainL 

aA 


H  0  implies  that  the  sum  of  the  “refined  residuals”  should  be  zero.  Solv¬ 


ing 


ainL 

as 


=  0  implies  that  the  vector  of  refined  residuals  should  be  orthogonal  to  the 


vectors  of  means.  It  follows  that  the  optima!  estimate  of  A  and  B  can  be  obtained  by 
the  “standard  least  squares”  inversion  with  the  censored  oata  all  replaced  by  their 
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conditional  expectations,  i.e. ,  the  “refined  observations".  Thus  a  can  be  solved  itera¬ 
tively  with  [6]  along  with  a  and  p  using  the  EM  algorithm.  In  the  non-censored  case, 
this  “MLE-CY"  code  gives  results  identical  to  those  derived  by  the  standard  least 
squares. 


If.3.3  A  General  m^-Yield  Regression  Routine  with  Uncertain  Data:  DWLSQ 

Even  the  19  Semipalatinsk  explosions  for  which  the  “exact"  yields  were  published 
would  inevitably  be  subject  to  many  sources  of  error.  The  Soviets  might  have  rounded 
8  of  the  announced  19  yields  to  the  nearest  5  kt  or  10  kt.  An  announced  yield  of  100 
kt  (e.g. ,  660320D)  could  mean  something  actually  measured  between  95  kt  and  104 
kt.  It  could  also  indicate  that  100  kt  was  the  designed  energy  release,  and  the  actual 
yield  was  somewhere  nearby.  Likewise,  the  "real  yield"  of  2  kt  {e.g. ,  720902M)  could 
be  something  between  1.5  kt  and  2.4  kt.  Below  100  kt,  the  rounding  errors  could 
overwhelm  the  presumed  standard  measurement  error  —  assuming  the  announced 
yields  are  not  otherwise  “fudged”. 

A  more  general  regression  routine  has  been  developed  to  take  the  rounding  and 
standard  errors  in  the  yields  into  account  (Jih,  1991).  For  each  (m^,  yield)  pair,  we 
use  a  random  number  generator  to  produce  a  perturbed  (m^,  yield)  pair  according  to 
their  uncertainty  distribution.  A  standard  least-squared  regression  is  then  performed  for 
each  data  set  of  19  perturbed  pseudo-observations.  The  procedure  is  repeated  for 
several  hundred  iterations,  and  all  the  resulting  calibration  curves  are  then  used  to 
infer  the  ensemble  behavior.  This  “doubly-weighted  least-squares  scheme"  is  an 
extension  to  the  "ordinary  weighted  least-squares”  in  which  only  errors  in  the 
would  be  used  to  adjust  the  inferred  parameters. 

The  “upper  95%  confidence  limit”  of  the  predicted  at  a  given  log(yield)  level 
(say,  Yq)  can  be  computed  as  follows; 
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itiij(max)  +  t(D.O.F.,0.975)[c^(mi,)  +  a^(regression)(-^  + 


(Yq  -  Y)^ 
KYi  -  Y)2 


)j0.5 


(7] 


where  N  =  number  of  data  points  used  in  the  regression,  D.O.F.  =  N-2.  o{mb)  =  the 
mean  S.E.  in  the  network  mt,  used  in  the  regression,  a(regression)  =  the  o  of  residu¬ 
als,  mb(max)  =  estimate  of  the  largest  possible  mean  at  the  given  log{yield)  level. 
Y  is  the  mean  log(yield)  used  in  the  regression,  and  t(D.O.F.,  0.975)  is  the  97.5  per¬ 
centile  of  Student's  t  distribution  at  “D.O.F.”  degrees  of  freedom.  The  “lower  95% 
confidence  limit”  can  be  computed  in  a  similar  way. 


ii.3.4  Expansion  of  Geotech’s  WWSSN  Database 

Our  database  of  station  values  based  on  the  short-period  vertical-component 
(SPZ)  WWSSN  [World  Wide  Standard  Seismograph  Network]  recordings  of  body 
waves  has  been  expanded  to  192^  events  from  a  variety  of  regions  including  the 
N.T.S.  (U.S.).  French  Sahara,  Azgir  (U.S.S.R.),  Urals  (U.S.S.R.),  Murzhik  (E.  Kazakh, 
U.S.S.R.),  Degelen  Mountain  (E.  Kazakh,  U.S.S.R.),  Balapan  (E.  Kazakh,  U.S.S.R.), 
Novaya  Zemlya  (U.S.S.R.),  Tuamoto  Islands  (France),  Rajasthan  (India),  and  Lop  Nor 
(Sinjiang,  China).  This  database  consists  of  515  usable  "a”  (/.e. ,  zero-crossing  to  first 
peak),  “b”  (/.e. ,  first  peak  to  first  trough),  and  “max”  (/.e. ,  max  peak-to-trough  or 
trough-to-peak  in  the  first  5  seconds)  event  phases.  Jih  et  al.  (1990b)  reported  evi¬ 
dences  which  indicate  that  the  GLM  network  values  interred  from  these  WWSSN 
recordings  are  better  than  many  other  magnitude  measurements. 


’Ill  events  were  measured  by  R  A  Wagner,  M  E  Marshall.  R  O  Ahner,  and  J  A  Burnetti  under  previous  contracts  dur¬ 
ing  the  past  decade  R.  A  Wagner  added  56  more  events  to  Geotech's  m^^  database  during  FY90'91  The  remaining  25  events 
were  adapted  from  the  data  that  DARPA  distributed  in  1988 
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11.3.5  Magnitude-Yield  Relationship  at  Semipalatinsk  Area 

A  systematic  comparative  analysis  of  the  magnitude-yield  relationship  at  three 
regions  (Balapan  [=  Shagan  River]  ,  Degelen,  and  Murzhik  [=  Konystan])  of  Semipala¬ 
tinsk,  Eastern  Kazakhstan,  U.S.S.R.  has  been  conducted  using  miscellaneous 
unclassified  magnitudes  as  well  as  the  recently  published  yields  of  96  Soviet  explo¬ 
sions.  Only  the  most  noteworthy  observations  are  summarized  here: 

•  1  MLE-CY  and  DWLSQ  vs.  Ericsson’s  code.  Including  the  censored  yields  in  the 

regression  does  generally  improve  the  accuracy  of  the  yield  estimates  slightly,  if 
the  magnitudes  are  “consistent".  The  “MLE-CY"  code  is  robust  in  detecting  the 
inherent  inconsistency  of  the  magnitudes  by  utilizing  the  censored  information. 

In  reality,  both  the  magnitude  and  the  yield  measurements  are  subject  to  error. 
Pending  the  determination  as  to  which  of  the  two  extreme  hypotheses,  namely 
a(mb)/a(Y)  =  0  and  0(mb)/a(Y)  =  oo,  is  closer  to  the  real  situation,  we  also  applied 
Ericsson’s  (1971)  curve-fitting  method^  to  the  regressions  of  non-censored  yields 
using  various  a(mfa)/a(Y)  ratios.  As  expected,  we  can  see  the  smooth  transition  of 
estimated  parameters  (i.e.,  the  slope  and  the  intercept)  as  a(m^)/a{Y)  varies. 
Thus  the  censored  cases  with  nontrivial  o(mb)/o{Y)  values  could  also  be  “interpo¬ 
lated”  accordingly  (Jih  et  al. ,  1990b).  Note  that  Ericsson’s  (1971)  method  allows 
different  variances  in  both  the  independent  and  the  dependent  variables  in  the 
regression.  However,  it  can  be  applied  to  the  non-censored  case  only.  "MLE-CY” 
and  Ericsson's  methods  represent  two  different  directions  in  extending  the  stan¬ 
dard  least  squares.  “DWLSQ”  is  even  more  flexible  than  Ericsson’s  code  in  that 
it  permits  the  errors  in  each  (m^,  yield)  pair  to  be  arbitrary. 

•  2  Rounding  Errors  vs.  Gaussian  Errors  in  the  Yields.  There  has  been  some 

concern  about  the  accuracy  and  precision  of  the  Soviet  published  yields.  It  turns 
out  that,  so  long  as  the  best  (such  as  g )  is  used,  the  uncertainty  factor  in 
the  predicted  yield  of  future  Semipalatinsk  events  is  not  very  sensitive  to  the 

‘^Codo  provided  R  M  Shurriwny  and  T  M  McEilrosh 
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postulated  uncertainty  (precision)  in  the  published  yields  of  the  19  events.  Also, 
the  yields  of  future  underground  explosions  in  the  Semipalatinsk  area  can  be 
estimated  seismically  with  a  capability  much  better  than  the  factor-of-2  uncertainty 
that  is  commonly  reported.  For  instance,  a  factor  of  1.5  could  be  a  reasonable 
uncertainty  estimate  at  around  the  50-kt  level  if  9  is  used  as  the  source  meas¬ 
ure.  At  yields  below  10  kt  small  variations  of  the  physical  environment  may  pro¬ 
duce  greater  uncertainty.  Therefore,  the  uncertainty  may  be  inherently  greater  at 
such  low  yield  level. 


Table  4.  Inferred  Calibration  Parameters  and  Associated  Uncertainty  Factors 

Uncertainty  in  yield 

Slope 

Intercept 

1 

Ikt 

_ 

lOkt 

50kt 

lOOkt 

150kt 

tmmam 

R.E.  +  20%  S.E. 

0.793±0.031 

4.30810.050 

4.14 

2.09 

2.55 

2.89 

R.E.  +  10%  S.E. 

0.804±0.022 

4.28810.038 

3.23 

1.67 

1.88 

2.04 

R.E.  +  5%  S.E. 

0.805±0.020 

4.28610.035 

3.11 

1.78 

1.53 

1.64 

1.77 

R.E.  +  2%  S.E. 

0.806±0.020 

4.28510.035 

3.33 

1.81 

1.53 

1 

1.69 

1.82 

R.E.  +  1%  S.E. 

0.805±0.020 

4.28710.035 

3.04 

1.80 

1.49 

1.60 

1.75 

R.E.  Only 

0.80710.019 

4.28410.034 

1.83 

1.54 

1.60 

1.72 

20%  S.E. 

0.79410.031 

4.30610.049 

4.08 

2.16 

1.88 

2.23 

2.50 

10%  S.E. 

0.80710.018 

4.28210.027 

2.83 

1.74 

1.60 

1.77 

1.91 

5%  S.E. 

0.81110.012 

4.27710.018 

2.41 

1.62 

1.51 

1.61 

1.71 

2%  S.E. 

0.81110.009 

4.27610.014 

2.33 

1.57 

1.48 

1.57 

1.64 

1%  S.E. 

0.81210.009 

4.27510.014 

2.32 

1.57 

1.49 

1.57 

1.64 

0.1%  S.E. 

0.81210.009 

4.275+0  014 

2.32 

1.57 

1.49 

1  57 

1.64 

R  t.;  Rounding  Error  ,  S.E.:  Standard  Error 


•  3  Yield  estimates  of  recent  Balapan  explosions.  Ringdal  and  Hokland  (1987) 
noted  that,  for  Balapan  explosions  detonated  after  1976,  the  largest  peak  of 
clustered  NORSAR  RMS  Lg  values  was  6.06.  It  corresponds  to  138  kt  on  the 
calibration  curve  derived  with  “DWLSQ”: 


RMS  Lg  =  4.54{±0.050)  +  0.71  (±0.028)  log(W)  [8] 

which  is  almost  identical  to  the  mean  yield  of  139±7  kt  computed  by  Sykes  and 
Ruggi  (1989)  using  the  9  largest  Balapan  shots  during  1976-1985. 

•  4  Depth-yield  scaling  at  Eastern  Kazakhstan.  The  Soviet-announced  burial 
depths  [DOB]  indicate  a  strong  tendency  to  detonate  Semipalatinsk  explosions  at 
a  “scale  depth”,  Ds,  {i.e. ,  the  DOB  scaled  to  a  yield  of  1  kt)  of  117  meters,  based 
on  the  relation 

DOB  (meters)  =  117  •  [W(kt)]°-25  [9] 

[9]  is  determined  using  18  Semipalatinsk  events  (4  Balapan,  6  Murzhik,  and  8 
Degelen)  of  known  yields  and  DOB.  The  yields  are  assumed  to  be  perturbed  by 
10%  standard  error,  and  the  DOB  are  subject  to  0.1%  error.  Deleting  the  8 
Degelen  events  gives  very  similar  result  with  smaller  scatter  and  slightly  larger 
Ds; 


DOB  (meters)  =  145  •  [W(kt)]°-2'‘  [10] 

For  explosions  from  Murzhik  and  Balapan  test  sites,  the  focal  depths  correlate 
with  the  announced  yields  even  better  than  m5(NEIS),  mb(EKA),  mf,{4  UK 
arrays),  and  log(/Wo,  4  UK  arrays)  {cf.  Table  5D  of  Jih  etal.,  1990b).  Both  [9] 
and  [10]  strongly  suggest  that  the  DOBs  of  historical  Semipalatinsk  explosions 
appeared  to  be  proportional  to  the  quartic  root  of  the  yields  instead  of  the  cubic 
root  as  observed  at  NTS.  Degelen  explosions  tend  to  be  underburied  except  the 
event  660507.  The  Soviets  seem  to  have  been  pushing  the  DOBs  to  the  shallow 
limit  (and  thereby  accepting  the  possible  containment  risks)  at  Degelen  Mountain 
during  the  period  1961-1972. 
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•  5  Test  site  bias.  There  have  been  several  studies  which  indicate  that  Nuttli’s 
(1987)  “Degelen  puzzle”  (Table  5)  could  be  invalid  because  of  the  relatively 
poorer  quality  mi,{\SC)  used  (jih  etai,  1990).  The  updated  bias  estimate 
between  Eastern  Kazakhstan  and  NTS  is  systematically  larger  than  0.35  if 
values  reported  by  TG,  S-cubed,  or  UK/AWE  are  used  (Table  6;  see  also 
Evernden  and  Marsh,  1987).  The  WWSSN  data  reveal  a  bias  of  0.13  m.u. 
(relative  to  NORSAR’s  RMS  Lg)  between  the  SW  and  NE  subregions  of  Balapan 
Test  Site,  which  confirms  what  Ringdal  and  Marshall  (1989)  found  with  ISC  and 
NORSAR  data.  Relative  to  m^,  -yield  relationship  does  appear  to  be  more 
transportable,  as  indicated  by  the  insignificant  difference  between  KTS  and  NTS 
calibration  curves  using  RMS  Lg  (Table  6). 


Table  5.  Nuttli’s  (1987)  Estimates  of  mt.  Bias  (Relative  to  rnij{Lg) 


Test  Sites 

Magnitude  Used 

10  kt 

100  kt 

150  kt 

Balapan  -  NTS 

rnb(lSC),  mi,{Lg) 

0.35 

0.35 

0.35 

Degelen  -  NTS 

r77b(ISC).  mt,{Lg) 

0.58 

0.58 

0.58 

-  82  - 


Table  6.  Updated  Estimates  of  Test  Site  Bias 

Magnitudes  Used 

KTS-NTS  Bias 

KTS 

NTS 

10  kt 

100  kt 

150  kt 

Marshall  et  al.  ^ 

Marshall  et  al.^ 

0.55 

0.48 

0.46 

Jih  and  Wagner^ 

Jih  and  Wagner^ 

0.47 

0.42 

0.41 

Jih  and  Wagner^ 

Murphy® 

0.36 

0.36 

0.36 

Murphy® 

Murphy® 

0.47 

0.41 

0.40 

Ringdal^ 

Patton® 

0.09 

0.05 

0.04 

Israelson® 

Patton® 

0.04 

0.06 

0.07 

1)  Combining  all  UK/AWE's  Galapan,  Oegelen,  and  Murzhik  values  as  listed  in  Vergino  (1989). 

2)  UK/AWE's  NTS  m^  values  as  distributed  in  1987. 

3)  m^j  =  0.81  log(W)  +  4.28  using  19  KTS  explosions  with  published  yields  (Jih  and  Wagner,  1991). 

4)  m^  =  0  86  log(W)  +  3  76  for  NTS  high-coupling  events  (Jih  and  Wagner,  1991). 

5)  mjj(S-cubed)  =  0.75  log(W)  +  4  45  based  on  the  network-averaged  spectra  (Murphy,  1990) 

6)  mjj(S-Cubed)  =  3.92  +  0.81  log(W)  for  NTS  high-coupling  events  (Murphy.  1981). 

7)  RMS  Lg  =  0.71  log(W)  +  4.54  with  NORSAR  RMS  furnished  by  Ringdal  (1990). 

8)  ma(Lg)  =  0.76  log(W)  +  4.40  with  LLN  data  (Patton,  1988). 

9)  RMS  Lg  =  0  78  log(W)  +  4  42  based  on  RMS  Lg  values  furnished  by  Israelson  (1991b)  (Jih,  1991). 

•  6  Cratering  to  non-cratering  correction.  The  new  calibration  curve  for  Balapan 

explosions  also  provi(jes  an  alternative  and  straightforward  approach  to  derive  the 
rriij  adjustment  converting  cratering  shots  to  contained  explosions  of  the  same 
yield  {cf.  Table  8A  of  Jih  etai,  1990b).  The  correction  derived  by  this  approach 
matches  that  by  other  studies  rather  well. 

•  7  Distinct  features  of  Degelen  and  Balapan  sites.  Degelen  Mountain  is  the  only 

test  site  that  has  a  decreasing  log(Pmax/^a)  with  increasing  yields 

{cf.  Table  8A  of  Jih  et  al. ,  1990b).  It  is  also  the  only  test  site  for  which  the  phase 
“a”  shows  the  smallest  scatter  around  the  calibration  curve,  as  compared  to  the 
phases  “b"  and  "max”  (cf.  Table  6C  of  Jih  et  al.,  1990b).  Both  the  mountainous 
topography  (which  causes  complex  pP  interference)  as  well  as  the  testing  prac¬ 
tice  (e.g. ,  the  abnormally  shallow  shot  depths  and  the  usage  of  tunnels)  could  be 
responsible.  At  Balapan,  the  phase  “b"  has  the  smallest  scatter  around  the 
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calibration  curve  {cf.  Table  5C  of  Jih  etai,  1990b).  These  observations  confirm 
the  conjecture  that  the  first  cycle  could  give  better  results  than  does  the  “max” 
phase  in  a  proper  environment. 

•  8  Benchmark  of  various  magnitudes.  For  Balapan  events,  RMS  Lg  reported  at 

NORSAR  (Ringdal  and  Marshall,  1989;  Ringdal  and  Hansen,  1989)  and  W29 
based  on  WWSSN  provide  the  smallest  scatter  around  the  calibration  curve.  In 
fact,  even  7772.2  would  seem  to  be  better  than  almost  all  other  unclassified  magni¬ 
tudes  based  on  the  teleseismic  P  waves  or  log('PJ  in  terms  of  yield  estimation 
as  well  as  the  scaling  against  Ringdal’s  RMS  Lg. 

•  9  Abnormal  events  detected  by  MLE-CY.  The  Balapan  cratering  event  650115 

(100-150  kt,  178  meters)  and  four  events  at  Degelen  Mountain  (641116,  194 
meters;  660629,  187  meters;  661019,  185  meters;  and  671017,  181  meters)  were 
rejected  by  the  “MLE-CY”  code  as  outliers.  All  4  of  these  Degelen  events  were 
said  to  have  yields  between  20  and  150  kt.  However,  a  fully  contained  explosion 
at  Balapan  or  Murzhik  regions  with  a  shallow  DOB  of  180  meters  or  so  would  be 
expected  to  have  a  yield  near  2  kt  {e.g. ,  Murzhik  event  720902)  rather  than  any 
value  between  20  and  1 50  kt.  Another  interesting  event  is  Degelen  event  660507 
which  has  an  announced  yield  of  4  kt  and  a  remarkably  deep  DOB  (as  compared 
to  balapan  and  Murzhik  explosions).  Perhaps  an  experiment  with  a  much  larger 
yield  was  planned  for  that  explosion. 


11.3.6  REPORTS,  PRESENTATIONS,  AND  PUBLICATIONS 
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Meeting,  Victoria,  British  Columbia.  Canada. 
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Inversion  in  Heterogeneous  Media  (July  31  -  August  3,  1989,  Boston,  MA.) 
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(1991)  A  refined  network  m^  determination  scheme  incorporating  near-source  effects, 
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CO.)  (Eds  J.  Lewkowicz  and  J.  McPhetres),  PL-TR-9 1-2208,  Phillips  Laboratory, 

Hanscom  Air  Force  base,  MA.  (ADA24132S) 

li.4.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  new  magnitude  determination  scheme  (Equation  [1])  significantly  reduces  the 
fluctuational  variation  across  the  recording  stations,  as  illustrated  in  this  study.  It  is 
shown  that,  by  applying  this  scheme  to  worldwide  explosions,  it  is  possible  to  have  a 
consistent  base  line  in  estimating  the  absolute  magnitudes  (which  is  crucial  in  estimat¬ 
ing  the  test  site  bias)  while  the  precision  in  the  resulting  network  values  can  be 
maintained  as  well  as  could  be  achieved  by  the  single-test-site  approach.  The  stan¬ 
dard  error  in  most  values  is  0.02  m.u.,  about  the  same  as  that  for  RMS  Lg 
inferred  from  in-country  regional  recordings  reported  by  Israelson  (1991a)  and  Hansen 
et  at.  (1990). 

The  most  detailed  description  of  the  wave  propagation  model  would  naturally  sug¬ 
gest  that  yet  another  term  could  be  added  into  the  Equation  [1]  to  count  for  the 
source-region  attenuation.  So  far  such  source  region  bias  in  has  always  been 
inferred  with  other  information  such  as  the  yields  (as  in  this  study).  Mg  (e.g. , 
Evernden  and  Marsh,  1987)  or  P„  velocity  (e.g. ,  Marshall  eta/..  1979)  eta.  The 
hypothesis  that  such  attenuation  differential  could  be  directly  discerned  from  mi,  alone 
by  further  improving  Equation  [1]  is  worth  testing. 

The  two  regression  routines  developed  in  this  study  {i.e. ,  “MLE-CY”  and 
“DWLSQ”)  represent  two  very  different  directions  in  extending  the  standard  least- 
squares  regression.  They  should  be  merged  together  for  a  more  general  and  robust 
tool. 

Digital  signals  recorded  on  in-country  seismic  instruments  at  regional  distance  will 
be  critical  for  monitoring  low  yield  explosions  below  10  kt.  Obviously  the  regression 
routines  developed  in  this  study  can  well  be  applied  to  other  source  measures  which 
are  based  on  regional  phases.  On  the  other  hand,  teleseismic  data  such  as  the 


WWSSN  data  used  in  this  study  still  carry  invaluable  information  that  are  worthy  of 
further  exploitation,  beyond  simply  calculating  the  “unified  yield”. 

Throughout  this  study,  our  emphasis  has  been  to  reveal  the  site-dependent 
characteristics  from  the  observations  exclusively  so  that  in  the  future  the  inferred 
results  can  be  critically  examined  and  compared  with  those  derived  by  other  means. 
We  suggest  that  the  follow-up  research  be  accompanied  by  well-constrained  forward 
modeling  studies  using  realistic  structures.  The  upgraded  LFD  code  which  incor¬ 
porates  boundary  conditions  for  topographical  free-surface  of  arbitrary  shape  (Jih 
et  al. ,  1988)  in  addition  to  the  “strain  filter”  and  “marching  grid”  features  as  outlined  in 
Jih  et  al.  (1989)  can  be  utilized  in  the  future  to  improve  our  understanding  of  the  fun¬ 
damental  issues  of  seismic  energy  partitioning  on  the  focal  sphere  as  well  as  their 
implications  for  yield  determination. 
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